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I.  Introduction 


Transataospherlc  vehicle  (National  Aero-Space  Plane)  technology  presents  a 
grand  challenge  of  the  rest  of  this  century  and  the  next  for  aerospace  science 
and  engineering  coaaunities  [l].  In  contrast  to  ballistic  re-entry  vehicles 
(space  shuttles)  piercing  through  the  ataosphere  nearly  vertically,  the  trans- 
ataospberic  vehicle  will  reaain  in  the  dense  layer  of  the  ataosphere  for  a 
good  portion  of  its  flight,  and  hence  the  air  near  the  aoving  vehicle  is  heat¬ 
ed  to  vory  high  teaperatures  (aerodynaaic  heating).  According  to  conputational 
fluid  dynaaic  siaulations  of  hypersonic  lifting  body,  the  skin  teaperature  can 
readily  reach  2000°F  and,  in  particular,  the  nose  cone  and  external  skin 
panels  near  the  jet  efflux  are  eotlaated  to  have  teaperatures  over  3500°F  [lj. 
At  these  hlgh-teaperature  hypersonic  flights,  acoustic  fatigue  becoaes  a 
severe  structural  problea  because  not  only  are  the  pressure  fluctuations  anti¬ 
cipated  in  the  range  of  160-180  dB,  but  also  the  skin  panels  can  vibrate 
chaotically  about  the  theraally  buckled  positions  whose  aaplitudes  increase  as 
the  square  root  of  teaperature  above  the  critical  buckling  teaperature  (oil¬ 
canning  effect)  [2].  Realizing  the  iaportance  of  elevated  teaperatures  on 
acoustic  fatigue,  work  has  begun  as  early  as  in  the  70 's  to  establish  fatigue 
failure  criteria  in  a  contained  tberaal-acoustlc  envlronaent  [3]. 

Strictly  speaking,  sonic  fatigue  is  a  fluid-structure  interaction  problea 
which  requires  slaultaneous  solution  of  the  Navier-Stokes  equations  for 
pressure  and  teaperature  together  with  the  structural  dynaaic  equations  for 
ensuing  stress/strain  distributions.  For  large-aaplitude  deflections,  the 
shape  of  structures  representing  the  fluid-structure  interface  is  not  known  a 
priori  ,  and  hence  aust  be  deternined  fron  slaultaneous  solution  of  the  fluid 
and  structural  equations.  However,  auch  progress  has  been  aade  in  the  past  by 
decoupling  the  structural  part  of  the  problea  froa  solving  the  Navier-Stokes 
equations  for  pressure  fluctuations.  Following  in  this  tradition,  we  shall 
in  this  report  investigate  dynaaics  of  a  piece  of  hypersonic  vehicle  structure 
aodelled  by  the  so-called  von  Karaan-Herrnann-Chu  plate  equations  subjected  to 
prescribed  pressure  fluctuations  and  teaperature  variations. 

Although  Bolotin  [4]  derived  the  basic  plate/shell  equations  for  large- 
aaplitude  deflection  subjected  to  teaperature  variation  in  the  late  60's  (Sec. 
II).  it  is  fair  to  say  the  role  of  theraal  effects  has  not  yet  been  elucidated 
in  its  full  generality.  This  is  perhaps  due  to  the  coaplexlty  of  nonlinear 
equations  and,  aoreover,  a  wide  variety  of  boundary  and  edge  conditions  that 
one  aay  lapose  on  the  plate.  To  exhibit  the  essential  physics,  wa  expand  the 
transverse  dlsplaceaent,  Airy's  stress  function,  and  teaperature  distribution 
in  trigonoaetric  functions,  and  thus  obtain  acdal  equations  for  both  siaply- 
supported  and  claaped  plates  (Sec.  III).  Note  that  the  original  plate  equa¬ 
tions  are  partial  differential  equations,  whereas  the  nodal  equations  are 
ordinary  differential  equations.  However,  the  price  pajd  for  thia  reduction 
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(Galerkin's  procedure)  is  lndeterminancy  because  the  system  of  nodal  equations 
is  not  closed  because  of  nonlinearity.  We  shall  therefore  truncate  the  system 
to  obtain  modal  equations  for  the  first  four  even  modes  in  Sec.  IV,  and  there¬ 
by  permitting  comparison  with  the  previous  formulations  of  Levy  [5]  and  Paul 
[6],  It  is  important  to  point  out  that  after  a  suitable  nondimenslonallzation 
the  modal  equations  of  simply-supported  and  clamped  plate  can  exhibit  a  simi¬ 
lar  form  for  the  three  temperature  terms.  The  flr9t  Is  global  thermal  expan¬ 
sion  by  uniform  temperature,  the  second  corresponds  to  local  thermal  expansion 
by  temperature  variation  over  the  plate,  and  the  third  term  represents  the 
thermal  moment  owing  to  temperature  gradient  through  the  plate  thickness. 

Even  for  the  low  order  modal  equations,  there  are  too  many  cubic  terras  for 
us  to  readily  ascertain  the  correctness  of  Galerkin's  procedure  carried  out  in 
Sec.  IV.  We  therefore  demonstrate  credibility  of  the  modal  equations  by  way 
of  constructing  the  Hamiltonian  which  embodies  the  kinetic  and  potential 
(strain)  energies  of  the  plate  being  conserved  in  the  absence  of  viscous  damp¬ 
ing  (Sec.  V).  To  proceed  further,  it  is  necessary  to  specify  the  temperature 
variation  over  the  plate  and  temperature  gradient  across  the  plate.  We  have 
adopted  in  Sec.  VI  very  simple,  but  nontrivial  profiles  for  the  temperature 
variation  and  gradient  to  expedite  the  subsequent  analysis. 

Much  insight  into  thermal  terms  can  be  gained  by  the  prototype  single-mode 
equation  shared  by  both  the  simply-supported  and  clamped  plates  (Sec.  VII). 
First  of  all,  the  uniform  temperature  and  temperature  variation  represent  a 
kind  of  thermal  stiffness,  but  they  add  negatively  to  the  structural  stiff¬ 
ness.  Hence,  the  combined  thermal-structural  stiffness  remains  positive  when 
the  sum  of  uniform  and  local  temperatures  is  less  than  the  critical  buckling 
temperature  (pre-buckllng) ,  whereas  It  becomes  negative  when  the  sum  exceeds 
the  critical  buckling  temperature  (post-buckling).  In  contrast,  the  tempera¬ 
ture  gradient  across  the  plate  gives  rise  to  thermal  moment,  hence  represent¬ 
ing  an  additional  loading,  as  already  observed  by  Boley  and  Weiner  [7].  There¬ 
fore,  this  together  with  external  pressure  forces  constitute  the  combined 
thermal-applied  forcing. 

For  the  acoustic  loading  it  i9  necessary  to  consider  stochastic  dynamic 
formulation  of  the  single-mode  equation  and  estimate  the  mean  square  response 
amplitude  subjected  to  Gaussian  random  excitations  [8].  Although  the  equi¬ 
valent  linearization  technique  [9]  has  proven  useful  for  nonlinear  structural 
dynamic  problems,  it  cannot  be  applied  directly  to  the  present  problem  because 
the  thermal  moment  appears  as  an  additional  time-independent  forcing.  This 
therefore  calls  for  an  extension  of  the  equivalent  linearization  to  nonzero- 
nean  Gaussian  excitations.  The  main  thrust  of  this  report  (Sec.  VIII)  Is  to 
show  the  growth  of  mean  square  displacement  as  the  plate  goes  through  thermal 
buckling,  the  competing  mechanical  and  thermal  loading,  and  the  effect  of 
various  thermal  terms  on  the  extreme-fiber  stress  and  strain  tensor 
components . 
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Most  detailed  information  hat  been  relegated  to  the  appendices,  for  the 
benefit  of  those  readers  who  demand  proofs.  In  any  event,  this  report  contains 
a  complete  Galer’;ln  formulation  of  simply-supported  and  clamped  plates, 
including  the  combined  stiffness  and  applied  forcing  terms. 


II.  Plate  Equations  for  Large-Amplitude  Displacement 

Let  us  begin  with  the  following  strain-stress  relations  for  a  plate, 
including  the  effect  of  thermal  expansion  «T 


cx  "  I  ''V  +  aT' 
*»  ”  i  <®y"  <"«>  *  “T' 
*xy’  «xr 


where  E  is  the  modulus  of  elasticity  and/*  Poisson's  ratio.  Besides,  T 
denotes  the  local  temperature  of  plate  with  the  thermal  expansion  coefficient 
«.  In  the  absence  of  «T,  Eq.  (1)  is  the  usual  linear  relationships  of  the 
strain  tensor  components  ex,  ey,  and  exy  and  stress  tensor  components  <rx,  «y. 
and  o_..  For  a  positive  a,  raising  T  would  simply  result  in  increased  strain, 
in  conformity  with  the  intuitive  notion  of  thermal  expansion.  By  solving  Eq. 
(1)  for  the  stress  tensor,  the  inverse  relation  is 

a  *  — ®T"t*x+  “  (1+m)«T], 

x  (1 -Ma)  x  y 

°v  *  7  E-a:  [gv+  ~  U+M)«T], 
y  (1-/**)  y 

axy“  Tifey  *xy*  (2) 

It  is  important  to  notice  the  negative  sign  for  «T  terms;  hence,  the  stress 
may  in  fact  decrease  as  T  is  raised. 

Following  Bolotin  [4),  we  decompose  T  Into 

T(x,y,z)  =  T(x.y)  +  z9(x,y),  (3) 


where  Tfx.yJ-h-1/  T(x,y,z)dz  is  the  temperature  averaged  over  the  plate 
~h/a 

thickness  h.  and  d(x.y)  is  the  temperature  gradient  across  h  (Fig.  1).  Note 
that  only  the  linear  tenperature  differential  in  z  is  Included  in  Eq.  (3); 
any  nonlinear  temperature  variations  are  ignored  according  to  the  thin  plate 
theory.  He  shall  first  outline  briefly  derivation  of  the  compatibility,  trans¬ 
verse  displacement  equation,  and  plate  edge  conditions. 

Compatibility  Equation 

For  the  mean  strain  tensor  *x,  *y  and  *xy  at  the  mid-plate;  l.e., 

,  h/a  _  h/a  „  h/a 

i«h  ‘/  c  dz,  t  «h  l/  t,Tdz,  and  *  -h  l/  *_,dz,  we  have  the  following 
*  -h/a  *  y  -b/a  y  y  -h/a  y 

expressions 


-4- 


(4) 


V  fc  ♦  *<§?>*• 

¥x,"  i<8»  ♦ls>  * 


deduced  fro*  geoaetric  considerations.  The  condition  of  strain  coapatlblllty 
Is  obtained  by  eliminating  u  and  v  fros  Eq.  (4)  through  cross  differentiation 


(5) 


Now,  by  averaging  Eq.  (1)  over  h  we  obtain  the  alternate  expressions  for  the 
Bean  strain  tensor 


*x  "  EH  ^Rx“  + 

*y  ■  E5  <V  "V  *  **• 

V  sf  V  <*> 


h/s  b/s  h/s 

Here,  N »/  o_dx,  N  »/  o„dz,  and  N  »/  o_vdz  are  forces  per  unit  length 
*  -h/s  *  y  -h/s  y  *y  -h/s  *y 

of  plate.  They  are  often  expressed  by  Airy's  stress  function  F 


N  ,  d*.F 
*  »y*' 


n  .  ai£. 

»  a** 


m  .  a** 

xy  Sx3y’ 


(7) 


which  autoaatically  satisfy  the  stress  equlllbrlua  at  the  aid-plate;  l.e., 
aNx/^x-*dNxy/by-0  and  bN^/ax+aNy/by-O.  Substituting  Eq.  (6-7)  into  Eq.  (5) 
yields  the  coapatlblllty  condition 


y‘p  .  eh.v*T  -  Eb((|^)*-  <^*>(^*>). 

where  V4-  V*V*  is  the  biharaonic  operator. 


(8) 


Transverse  Dlsplaceaent  Equation 

Next,  we  consider  the  balance  equation  for  the  shear  forces  and  noraal 
loading  q,  on  one  hand,  and  for  the  shear  forces  and  bending  aoaents  ax  and 
By,  and  twisting  aoaents  ■Xy"AyX>  on  the  other  hand.  Upon  elialnating  the 
shear  forces  froa  such  force  balance  equations,  we  obtain 


8 

- *  + 


Ss8 


+  if  ♦  (i  +  2n  fe-  ♦  a  -o 
dyS  *  515y  "x  y  *y  SlBy  qn  u* 


8y* 


(») 


relating  the  aoaents  with  the  noraal  loading  qL.  Me  first  express  the  aoaents 
in  teras  of  w 
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ft 

ft 


X 

y 


-D[  &  *  M 


3x* 

3aw 

3V4 


-D[  *  M 


3aw  + 

3y* 

aaw  + 

3x* 


«{ l+/*)0] , 
a(l+/i)9]  , 


*xy-  "D(l  Effigy’ 

where  D=Eh3/12(l-/i3 )  is  the  flexural  rigidity.  We  then  include  in  <J„the 
Inertial  force,  viscous  dawping,  and  external  pressure  p;  i.e., 

qn  =  "ph  '  ph*  It"  +  p* 


(10) 


(ID 


where  p  is  the  cross-sectional  mass  density  and  £  the  damping  coefficient. 

Note  that  ph£(3w/3t)  has  been  introduced  as  a  symbolic  representation  for 
viscous  damping.  Perhaps,  a  aore  practical  damping  model  would  be  £D74(3w/3t) 
of  Maekawa  ilOi-  Upon  inserting  Eqs.  (7),  (10)  and  (11)  into  Eq.  (9),  we 
obtain  the  equation  for  transverse  displacement  consistent  with  the  von  Karaan 
type  of  large-aaplitude  deflection 


P^-~  +  />h£|^r-  -  p  +  W4w  +  «(1^)W*£  «  .  (12) 

3t?  3x*  3y*  3y*  3x*  3j®y3*3y 


The  pair  of  Eqs.  (8)  and  (12)  is  that  given  by  Eqs.  (4.131)  and  (4.132)  (in 
which  kx  and  ky  are  set  to  zero)  of  Bolotin  [4],  and  also  agrees  in  fora  with 
Eqs.  (13.7.1)  and  (13.11.3)  of  Boley  and  Weiner  [7]. 


Plate  Edge  Conditions 

It  aust  be  pointed  out  that  the  conpatibility  is  a  statement  about  force 
balances  at  the  aid-plate.  Hence,  it  iaplies  certain  constraints  on  u  and  v. 
which  dictate  the  aoveaent  of  plate  edges.  Let  us  assume  that  the  solution  of 
Eq.  (8)  is  made  up  of  the  particular  solution  and  a  hoaog  neous  solution  of 
V*v=0  which  is  given  by 


F 


h 


P  x» 

~T~  - 


Vy’ 


(13) 


where  integration  contents  Px,  Py.  and  P^  represent  the  aeabrane  stresses. 

As  shown  in  Sec.  II  of  Ref.  [11],  the  immovable  edge  conditions  of  zero  in- 
planc  displacement  are 

5§y  -°  and  JSlo(lx>dxdy  “°  at  x“°*  a* 

^y  "°  Md  Jj{My>dxdV  *°  at  y-°*  b*  (14) 

Here,  the  integrals  of  displacements  along  the  plate  edge  are  suppressed  in  an 
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average  sense.  By  coabining  Eqs.  (4)  and  (6),  we  find  appropriate  expressions 
for  the  Integrands  in  teras  of  F  and  w 


—  -MS*- 

ty  -<*g|  *aT  ~  i<Sy>a-  (15> 


as  given  by  Bolotin's  Eq.  (4.140)  in  Ref  [4J.  On  the  other  hand,  the  novable 
edge  conditions 


and  u=  constant  at  x=*0,  a, 
and  Vs  constant  at  y=0,  b. 


(16) 


perait  free  aoveaent  of  edges  with  zero  inplane  stress.  Here,  again,  the 
vanishing  of  inplane  aeabrane  stress  is  iaposed  by  the  integral  constraint. 
We  shall  be  concerned  in  this  report  with  the  iaaovable  edge  conditions,  and 
thereby  exhibit  theraal  buckling  under  general  temperature  distributions. 


Stress  Tensor 

Rather  than  displaceaents,  the  stress  and  strain  are  the  aore  physically 
relevant  quantifiers  in  structural  analysis.  We  shall,  therefore,  present 
here  the  stress  coaponents  expressed  in  F  and  w,  froa  which  the  strain  compo¬ 
nents  can  be  recovered  by  Eq.  (1).  Let  us  begin  with  the  following  expressions 


Bu  _ 

_  daw 

X 

ox 

3x*  * 

3v 

-  3aw 

V 

5y  ‘ 

ay2’ 

lf3u 

.  avt 

xy" 

5l5y 

+  5xJ 

_  3aw 

z  53$y’ 


(17) 


computed  by  assuming  that  surfaces  which  are  parallel  and  noraal  to  the  aid- 
plate  remain  so  after  heating.  (Eq.  (17)  is  identical  to  Eq.  (12.2.1)  of 
Weiner  and  Boley  [7]  with  the  replaceaent  7Xy"2cxy.)  Let  us  insert  Eq.  (17) 
into  Eq.  (3)  to  obtain 


o 


xy- 


-(l^)aTj.. 

-(i+m)«t}, 


(18) 
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When  the  stress  tensor  is  averaged  over  h,  one  finds  that  only  the  syaaetric 
teras  survive  and  the  odd  (unsyaaetric)  teras  drop  out.  Hence,  we  obtain  froa 
Eq.  (18) 

5y  *  /Jtx  "  "y  * 

%  ♦  B  ■  *4?*  v  «»» 

Now.  substituting  Eq.  (19)  back  into  Eq.  (18)  yields  the  stress  tensor 
expressed  in  F  and  w  (with  the  use  of  Eq.  (7)) 


3aF 

Ez  f3aw 

,3aW) 

EotzO 

°x  i 

(l-^^dx* 

3y»J 

TOT)’ 

3aF 

Ez  paw  , 

EazO 

°y  K 

3x* 

(l-pa)Uy* 

3^-1 

1  3aP 

Ez  3aw 

axy*~ 

H  3x3y 

(i+M)  5x5y 

» 

after  circuitous  substitutions. 

In  retropseet.  the  assospticn  that  the  derivatives  3u/3x.  3v/3y,  3aw/3 x*. 
....  are  constant  across  the  plate  was  essential  in  arriving  at  Eq.  (19). 
Obviously,  this  cannot  be  true  in  coaposite  plates,  as  pointed  out  to  ae  by 
Steve  Whitehouse.  Furtheraore.  it  is  because  of  this  assumption  that  the 
nonlinear  strain  tensor  (i.e.,  the  teras  (2/2)(3w/&x)a,  (i/2)(3w/?y)a,  and 
(l/2)(8w/3x)(3w/3y)  being  incorporated  into  the  right-hand  sides  of  ex,  ey, 
and  «x y,  respectively)  also  gives  rise  to  the  suae  Eq.  (20).  as  pointed  out 
to  ae  by  Jay  Lavraea. 

Note  that  Eq.  (20)  is  identical  to  Bqs.  (29-31)  of  Choi  and  Vaicaitis  [12] 
when  Nx  and  Ny  are  set  to  zero  in  their  equations.  For  later  reference  we 
point  out  that  in  Eq.  (20)  the  first  tera  is  the  arabrane  stress,  the  second 
tera  the  bending  stress,  and  the  third  tera  represents  the  theraal  stress 
Induced  by  teaperature  gradient  across  h. 

Noraallzed  Coordinates 

It  is  convenient  to  scale  x  and  y  by  the  respective  sides  a  and  b  of  the 
plate  (Fig  1).  Introducing  x-x/a  and  y**y/b  and  subsequently  dropping  the 
karats,  we  have  the  pair  of  compatibility  and  displaceaent  equations  (/f«b/a) 


(21) 


(22) 


P  ♦ 


W' 


4 34W 

Bx4 


♦  a^H.1 

By4  J 


£ 


■afi  af? 

3\*  By^ 


+ 


a!s  a!*  _  2^h_  aio , 

a/1  Bj^  3«y  B^y* ' 


together  with  the  immovable  edge  conditions 


/wM$  -  '^3  *  -  fr<ls>Vdy  ■<>• 

-  <■^9  *  •**  -  I(l5>1}d’'dy  -°-  <23) 


to  be  considered  in  the  present  investigation.  And,  the  stress  tensor  now  has 
the  fora 


,  _i_  Ez  ftfs3aw  +  3aW| 

51  h*By*  l^(lV)1  By*J 

,2La!i _ is _ i^a  ♦ 

y  bba  Bx*  b^l-j^By* 


BszO 

mrj- 


Egzg 

TT^* 


(24) 


Because  of  the  coordinate  normal ization,  it  should  be  reaeabered  henceforth 
that  x  and  y  extend  over  (0,1). 


Boundary  Conditions 

Besides  the  edge  conrtrsints,  we  further  assure  that  the  plate  edges 
undergo  no  transverse  displaceaent 

w(0,y)«  w(l,y)»  w(x,0)»  w(x,l)«  0.  (25) 


For  a  simply-supported  plate  the  tangential  coaponents  of  the  bending  aoaent 
being  zero  implies  (Eq.  (12.4.2)  of  Boley  and  Weiner  (7]) 

♦  «ba(l+p)0  -  0  at  x-0,  1, 

Bx2 

£ 

♦  «ba(l+/0$  -  0  at  y-0.  1.  (26) 

By* 

Since  Eq.  (26)  is  the  inhomogeneous  boundary  conditions,  the  usual  sine 
expansion 

"U’V)  -  CA,  *«n  Vx)Vy>’  (27) 

where  $||(x)asin(Birx),  cannot  be  used  unless  6-0  around  the  plate  edge.  Then, 
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Eq.  (26)  degenerates  to  the  homogeneous  boundary  conditions 


*  0  at  x-0,  1  and  2_S?  «  0  at  y»0,  1.  (28) 

Bx2  By* 

On  the  other  hand,  for  a  clasped  plate  we  have  (Eq.  (12.3.1)  of  Boley  and 
Weiner  [7]) 

*  0  at  x-0,  1  and  *  0  at  y-0,  1,  (29) 

as  given  by  Eq.  (12.3.1)  of  Boley  and  Weiner  [7],  which  are  independent  of  the 
temperature  gradient.  We  say  therefore  continue  to  use  the  trigonometric 
function  expansion 

w(x’y>  “  CoCo"mnVx)*n<y>’  <30> 

where  ♦||(x)"COs(m+l)nx-cos(»-l)iTX,  as  has  already  been  applied  to  the  thermal 
buckling  [6]  and  sonic  fatigue  [20]  problems. 
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III.  The  Derivation  of  Modal  Equations 

Using  expansion  Sqs.  (2.27)  and  (2.30),  one  can  reduce  the  set  of  partial 
differential  Sqs.  (2.21-2.22)  to  ordinary  differential  equations  for 
However,  owing  to  the  presence  of  Airy's  stress  function,  this  reduction  would 
lead  to  an  Infinite  set  of  aodal  equations.  It  is,  therefore,  necessary  to 
truncate  the  aodal  equations  for  practical  coaputation  —  an  ad  hoc  procedure 
not  dictated  by  the  problea.  For  a  siaply-supported  plate  the  introduction  of 
Eq.  (2.27)  into  Eqs.  (2.21-2.23)  gives  rise  to  the  aodal  equations  with 
diagonal  aass,  daaping,  and  stiffness  aatrices.  On  the  other  hand,  Eq.  (2.30) 
will  result  in  the  aass,  daaping,  and  stiffness  aatrices  which  are  nondiagonal 
in  the  claaped-plate  case.  This  is  because  ♦  „  are  not  orthogonal  and  the  fora 
of  $a  changes  after  the  second  and  fourth-order  differentiations.  The  latter 
siaply  restates  that  Y^^fx^ty)  is  not  an  engenfunction  of  V4Y-XY-0.  We 
shall  first  present  in  Sec.  Ilia  the  aodal  equations  for  a  siaply-supported 
plate,  prior  to  the  aore  coapllcated  derive-  tlon  for  the  claaped  plate  in 
Sec.  Illb. 


A.  Siaply-Supported  Plate 

Expanding  the  particular  solution  of  Eq.  (2.8)  in  cosines,  we  write  in 
view  of  Eq.  (2.13) 


F  * 


pxbV 

5 — 


Py*8*8 

— 


+  Eh  if*  if*  F  cosprrx  cosqTty, 
p*o  q»o 


(1) 


where  Px  and  Py  are  integration  constants  for  the  hoaogeneous  biharaonic 
equation.  Here,  we  have  invoked  Pxy-0  in  anticipation  of  3aF/3X3y-0,  as 
required  by  the  iaaovable  edge  condition.  Also,  a  siailar  expansion  in 
cosines  is  assuaed  for  T  [13] 


T  *  t  +  if*  if*  t  cosprrx  cosqny,  (2) 

°  p*o  q**o 

where  t  is  constant  teaperature  over  the  plate.  Note  that  in  Eqs.  (1-2)  the 
tera  for  p-q-O  is  excluded  from  the  double  sub.  Instead  of  Eq.  (2.27),  let  us 
redefine  the  expansion  by 

wU.Y)  -  Cofi-a  "an  Vx)Vy)  *  (3) 

where  sin(anx)  are  the  orthonoraal  eigenfunctions. 

The  coBpatlbility  With  the  use  of  Eq.  (3)  the  right-hand  side  of  Eq.  (2.21) 
can  be  expanded  in  a  cosine  series 

I?*  SD  cosprrx  cosqny ,  (4) 

9X3y  3x*  3y*  p=o  q-o  pq 

where  8^  is  given  by  Eq.  (A7)  of  Appendix  A.  Hence,  after  substituting  Eqs. 
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(1-2)  into  Eq.  (2.21),  we  collect  the  following  coefficients  with  the  use  of 
Bq.  (4) 

F  - - -  .  (5) 

M  «•<*#•♦*)  (fe +  <?/*)* 

As  pointed  out  in  Appendix  A,  g__  consists  of  the  nine  subs  Bt  -  B  of  Levy 

pq  l  9 

£5J. 

The  displacement  In  parallel  to  Eq.  (3),  the  0  and  p  are  also  expanded  by  i>m 


For  the  derivation  of  nodal  equations,  it  is  convenient  to  put  Eq.  (2.22)  in  a 
symbolic  form 

Ri  Ra  +  **  +  °*  (7) 

where 

V  *  -w|f-  - 

R  «  Db"4f^4a^J  +  BP*-^-*--  ♦  ^“1  . 

*  *•  Bx4  3x*By*  By4  ' 

R,“  «(l+M)Db“*{p*^-|  +  , 

r  -  -  ^b-4^  aie  +  s£f  _  2  a 1*.  aiL.1 . 

4  '•Bx*  By*  By*  Bx*  BxBy  BxBy' 

We  then  Introduce  Eqs.  (1),  (3),  and  (6)  Into  Eq.  (7),  and  sort  out  the  com¬ 
ponents  for  w„„,  as  in  the  Gelerkin  procedure.  Because  of  the  orthogonality 

i  i  rs 

w*jdx“aj>  onP  can  write  down  at  once 

/J/J  RA(x)Vy)dxdy  m  +  Mzr1 '  Prs-  (8) 

/J/J  Rafr(x)fa(y)dxdy  »  ~(P*r*+  »a)a"r8.  (®) 

/J/J  R,fr(x)fs(y)dxdy  -  -  **)«„.  (10) 

On  the  other  hand,  the  treatment  for  R4  is  complicated  because  of  the  product 
term?;  in  F  and  w.  After  some  algebra,  we  obtain 

flJl  *,♦,<«>♦.<»>**»-  -  *„(»„•*„)•  <“) 

as  shown  by  Eq.  (B5)  of  Appendix  B.  Here,  we  have  indicated  the  explicit 
dependence  of  fr#  on  and  F^.  Note  that  Eq.  (11)  is  still  indeterminate 
owing  to  the  presence  of  constants  Px  and  Py ,  which  we  shall  evaluate  under 
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the  iaaovable  edge  condition  of  Eq.  (2.23).  Reaeabering  that  p«q»0  Is  excluded 
froa  the  suns  in  Eqs.  (1-2),  we  find  that 


v  -  *  <*<«!?>*»]  - 

V  -  {-^i('^2<<(ls>‘>>  *  «<!?>*»]  -  »-c)«t0}.  »2) 

where  «f  »-/*/*  fdxdy  is  the  average  of  f  over  the  plate.  After  evaluating 

«(|j)a>>  and  «(|^)a».  and  inserting  Px  and  Py  back  into  Eq.  (11),  we  write 
the  final  expression  in  the  following  fora 


/J/i  R*Vx)Vy)dxdy  *  1i  *  **  +  *a  +  V 


where 


T^Ehat 


1  ba(l-M) 
_  n^aEh 


T~  If*?  +  f  )"rs • 


V- 


k*  ~r«'  un 


(^*pa+qa)‘ 


I,-  .V„.  pST  -*.£.>  .  .‘(pT  «*.*  .  ,  V  n*«^)W.. 

*  2b4(l-i*ta)t  a.n»i  *°  a,n«i  1,11  a.n-t  *"  a,n«i  ■“  J  r* 


«  (w  _ ES _ > 

r*  «!*//»)* 


According  to  Eq.  (r- ),  by  setting  the  sua  of  Eqs.  (8-10)  and  (13)  to  zero  we 
obtain  the  aodal  equation  for  wrg  of  a  siaply-supported  plate.  Although  it  is 
possible  to  consolidate  Eq.  (13)  into  a  aore  coapact  fora,  we  prefer  to  leave 
it  ia  the  present  fora  involving  and  $r#,  for  readability  is  aore  iapor- 
tant  than  coapactness  wh^n  one  atteapts  to  enuaerate  the  aodal  equations. 


B.  Claaped  Plate 

Although  it  is  aost  desired  to  expand  w(x,y)  by  the  orthonoraal  eigen¬ 
functions  of  V4Y-XY*0  for  the  claaped-plate  boundaries,  the  use  of  such  eigen¬ 
functions  is  indeed  intractable  for  the  derivation  of  aodal  equations,  siaiiar 
to  what  we  have  done  in  Sec.  Ilia.  Therefore,  we  shall  be  content  here  with 
the  use  of  which  nay  be  expressed  alternately  by  ♦J>(x)*2S1|(x) ,  where 
SB(x)>sin(BTX)slmrx.  Since  lm  are  independent,  we  construct  the  orthonoraal 
functions  qB  froa  SB  by  the  Graa-Schaidt  procedure  [14].  Since  S^x).  S^tx), 
...  are  even  functions  of  x  and  SB(x),  Six),...  are  odd  functions,  one  finds 
that  the  orthonoraal ized <pB(x)  also  split  into  the  even-shaped  coaponents 

«P1(x)v€73si(x),  ^i(x)vi4/Ss#(x)vS7lftS1(x),  ... 
and  the  odd-shaped  coaponents 
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<Pa(x)*2Sa(x),  <p4  { x )  =v1«>/3S4  ( X )  +v*73Sz  ( X ) ,  ... 

As  att  illustration,  we  have  compared  in  Fie.  2  the  primitive  and  #3  with 
the  orthonormal ized  <pl  and  <?>3 .  Summing  up  the  even  and  odd  components,  we 
have  the  following  orthonormal  bases 


Vx)  *  amiSi(x)' 


where 


V/S73  0  o 

0 

0 

0 

0  2  0 

0 

0 

0 

X/S715  0 

0 

0 

0 

0  v^73  0 

vWT 

0 

0 

V§7W  0  3^6/35 

0 

v407T 

0 

0  0 

\/®/3 

0 

V* 

is  the  lower  triangular  matrix.  Hence,  in  contrast  to  Cq.  (2.30),  the 
orthonormal  expansion  for  a  clamped  plate  is  given  by 


»<»•»)  * 


Note  that  the  second  equality  gives  the  working  definition  for  modal 
expansion. 

The  compatibility  Substituting  Eqs.  (1-2)  and  (15)  into  Eq.  (2.21),  we  obtain 
by  collecting  the  coefficients  for  p  and  q 

«batM  P™ 

B  - - SSL -  + - M -  ( 16) 

Pq  «•<***♦*)  (P*)»  -  t/fi)* 

wt.ere  P__  is  given  by  Eq.  (C12)  of  Appendix  C.  Note  that  P__  is  considerably 
Ph  Pm 

more  complicated  than  8__  (Appendix  A)  due  to,  in  part,  the  use  of  orthonormal 

Ph 

bases.  However,  it  presents  no  serious  handicap  in  that  both  £__  and  J>  are 

Ph  P*i 

too  unwieldy  for  hand  enuweration  and  hence *  in  any  event,  will  have  to  he 
enumerated  by  symbolic  manipulations  on  a  computer  (see.  Appendix  D). 

To  establish  a  contact  with  Paul's  tabulated  results  (Appendix  B  of  Ref. 
[6]),  we  must  remember  that  his  analysis  was  based  on  the  nonorthogonal 
expansion  of  Bq.  (2.30).  However,  note  that  Eq.  (15)  does  reduce  to  Eq.  (2.30) 
under  ■jj-4  For  instance,  of  the  30  quadratic  terms  listed  in  Eq.  (Dl)  of 
Appendix  D,  there  are  10  terms  which  involve  only  the  diagonal  a^j.  Hence, 
with  a(l,l)«a(3,3)*l  the  10  terms  simplify  to  give 
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v48w.  -32w.  +160W*  -240w„« W„-80w*„ . 

13  31  13  33  31  31  33  33 


(17) 


On  the  other  hand,  Paul  presented  the  following  matrix  (the  tabulation  for 
P-G  and  Q-2  in  Appendix  B  of  Ref.  [6]) 


‘  32 

18 

-32 

-16' 

a 

-64 

16 

64 

-16 

-32 

-16 

160 

30 

.  34 

-16 

-320 

80. 

to  form  F„  by  KT5^.W,  where  W  la  the  column  vector  (w>4  ,w„,_ )  and  T  is 

the  transpose.  It  is  easily  checked  that  W  ft  W  is  identical  to  Eq.  (17), 
Furthermore,  other  components  of  F_,  can  also  he  checked  against  Paul's  tabu- 
lated  results  in  a  completely  analogous  manner. 

The  displacement  In  anclogy  to  Eq.  (6),  we  let 

and  carry  out  the  Galerkin  procedure.  In  view  of  /*<p j<P  jdx=6  j ’  we  have  by 
inspection 

/J/J  R1*r(x)*s<y>d*dy  *  pi^T  +  ’  prs‘  (19) 


For  the  biharmonic  term  we  obtain  after  some  algebra 
/£/£  *a9r(x)<Ps(y)dxdy  *  ^jr  V 


where  kt  is  given  by  Eq.  (E5)  of  Appendix  E.  Similary,  we  have  for  the 
temperature-gradient  term 

/£/£  R,er(x)<pa(y)dxdy«  -  (21) 

where 


Here,  the  summation  notation  Zjf (I)*f (l)+f (-1)  of  Maekawa  [10]  was  used  to 
consolidate  8  terms  into  the  2  terms  of  %  .  Finally,  the  Integral  involving 
R4  is  given  by  Eq.  (F6)  of  Appendix  F 
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*4*r(x*s<y)dxdys*  -  ^~(Pxb%+PyaaA9)  -  yrs<wBn*Ppq>-  (22> 


To  render  Eq.  (£2)  completely  determinable,  it  is  necessary  to  evaluate  Px  and 

Py  under  the  imaovable  edge  conditions.  Using  the  «{|^)*»  and  «  <#>  a» 
evaluated  in  Appendix  G,  we  again  rearrange  Eq.  (22)  into  linear  and  cubic 
contributions 


where 


f'0fo  R4<?r(x)^s(y)dxdy  "  Jt  +  Js  +  Js  +  V 
TT^Ehat.  „ 

J  -  -  - - —  +A. )  # 

4  b*(l~ii)  3  8 


J_=  - 


IT  (l-/i) 
naga«Eh 


ka  rs'  mn 


(0*P^ *9* ) 


V  M-,4»  *  ("V  "'“V4.}- 

J  -  -  s!£S!»  y  („  t  - l!fia - ) . 

♦  4b4  *rs'  an’  {^pa+ 


The  aodal  equations  for  a  clamped  plate  are  obtained  by  assembling  Eqs. 

(19-21)  and  (23)  according  to  Eq.  (7). 

For  completeness ,  attempts  were  made  in  Appendix  I  to  compare  ?rg  with  the 
corresponding  formulas  derived  by  Paul  [6].  However,  we  can  provide  only  a 
partial  comparison  because  the  orthogonality  of  Eq.  (14),  which  plays  an 
important  role  in  our  formulation,  is  completely  absent  in  Paul's  formulation. 


IV.  Lowest-Order  Nodal  Equations  of  Even  Node  Shape 

The  conplexity  of  Eqs.  (3.13)  and  (3.23)  prevents  us  froa  seeing  the 
overall  structure  of  nodal  equations.  Even  aore,  It  is  very  difficult  to 
assess  the  role  played  by  each  tera  and  thereby  deaonstrate  any  inter¬ 
relations  aaong  the  terns  of  nodal  equations.  To  this  end,  we  shall  in  this 
section  enuaerate  the  lowest-order  nodal  equations  which  include  only  the  four 
nodes  w,  .  w„„,  w„4  and  w„.  For  a  unified  representation,  however,  it  is 
necessary  to  put  the  nodal  equations  in  dinensionless  forn.  As  in  Ref.  [15], 
we  first  choose  the  length  scale  h,  the  tine  scale  7«[phb4/*4D]1/,a,  and  the 
force  scale  (ph?/7*).  We  then  introduce  the  tenperature  scale  T*,  which  will 
be  defined  later  for  the  sinply-supported  and  clanped  plates.  Using  these 
scales,  we  forn  the  dinensionless  aechanical  variables 


T  “t/7>  Prs  -<T*/Ph*)pr8.  Wrs  -wr8/h. 
and  thernal  variables 


T0  =t„/T», 


T  -t  /T*,  e  -he  /T*. 
pq  pq  rs  *rs' 


(1) 

(2) 


Because  of  the  factor  h,  the  0r#  now  represents  tenperature  differential 
rather  than  gradient  across  the  plate. 


A.  Siaply-Supported  Plate 

According  to  Eq.  (3.7),  the  nodal  equation  follows  by  equating  the  sun  of 

Eqs.  (3.8-3.10)  and  (3.13)  to  zero.  For  the  sinply-supported  plate,  the 

appropriate  T*  is  the  critical  buckling  tenperature  T*^t*i£  (/^+l)/13«l^  (1+p) . 

s 

at  which  the  global  thernal  expansion  cancels  out  the  aechanical  stiffness  for 
wii  Then,  the  dinensionless  nodal  equations  becones 


BaW 


dr 


rs  +  ^  a*r. 


Prs  +  ^^"rs 


(SSP-1)  (SSP-2)  (SSP-3)  (SSP-4) 

+6^a{r2(^aZ°°  *awjjL  +  plT  naW?  )+  ss(/i2f0  naN?  +  /Talf  naW*  ))w 
n,n=i  n,n«i  n,n-i  n,n-i  J 


(SSP-5a) 


-3^(1-Ma)f  (w  ,  - — £9- - - 

rs  ***  (^p8+  q*//>)a 

(SSP-5b) 


). 


umQeZni  ,n„m,  — ^3_) 

(SSP-6)  (SSP-7) 


-  er,  -  o. 

(SSP-8) 


(3) 


-17- 


We  may  interpret  the  terms  labeled  by  SSP  (simply-supported  plate)  1  -  8  as 
follows:  SSP-1  Is  the  inertial  tern,  SSP-2  represents  viscous  damping,  SSP--3 
is  the  external  forcing,  SSP-4  is  the  usual  stiffness  term,  and  terms  SSP-5a 
and  5b  represent  the  cubic  nonlinearity.  Note  that  SSP-5a  is  contribution 
from  the  iaaovable  edge  conditions,  whereas  SSP-5b  is  derived  froa  the  product 
teras  of  w  and  F.  Thermal  effects  are  embodied  by  the  last  three  terms.  That 
is,  SSP-6  is  the  global  thermal  expansion  owing  to  uniform  temperature,  SSP-7 
is  the  local  thermal  expansion  by  temperature  variation  over  the  plate,  and 
SSP-8  represents  the  thermal  moment  induced  by  temperature  gradient  across  the 
plate . 

Using  the  S__  and  Jt  tabulated  in  Appendices  A  and  B,  respectively,  Bq. 

Pm  rs 

(3)  for  r  and  s  «1  or  3  yields 


•  •  • 

W  +  r*W  -  P  +  a„  W  „  +  (SSP-5)  -  |^+l)TbM« 

rs  M  rs  rs  rs  rs  '  'P8  'r  9  o  rs 


o  rs  rs 


_  f  (..  Lml-  )  -  i^ilib  ©  -  0, 

4  rs'pSpa.j.qa  ’  12  rs°rs  ' 

where  the  overhead  dot  denotes  d/dr.  No  sum  is  implied  in  Eq.  (4)  by  the 
repeated  indices,  and 

ait“(0*+1)2’  aJl-(9)ra+l)a,  a^-SK^+l)*, 

b b3l-(9^+l),  baa-8(**+l). 

The  f  given  by  Eq.  (B6)  of  Appendix  B  are 


(4) 


=  hiW«+  Vx.+  h.*,i  +  VW 

*  Via*  Vxa*  he*ai+  V.,’ 

***{f^f)  *  h»Wii  +  Va.+  V*.i  +  V.*’ 

*aa^a ^  h7w1#+  b8W#1+  bi0W#J, 


(5) 

(6) 


(7) 


where 


V"2(Toa+Tao /fi*  }  *  ba-2T0a-4TM/{/?*rl)-2T04rTa4/(^+4)  , 

h,  -2Tao  ff  -4Taa /  (p*  +1 )  -2T40  /**  *T4a/  ( ♦  1 ) . 
h4=-9Ta4/(0*+4)  -9T4a/(^*+l),  18(Ta0/^*+T0#/9), 

h||-25Ta4/(8ar4)-25T4a/(4^+i)+18Taa/(^+l)+4T44/(^+1). 
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h7*18Tao/^*":l8T4o/^a"36Ta«/^+®)+9T*t/(4^+9)  •  *»,— W(Toa+Tto/9^,)  * 

h»“18Toa"18To4"36Taa/(9^8+1)+9Ta*/(9^a+4) *  hio“"18(Toa+Tao/^*>  *  (8) 

Finally,  the  components  of  SSP-S  are 


(JSP-5)^- 

(SSP-S),,- 

♦•.a8, *“,.",.’5. 

(SSP-S),,- 

(ssp-s),,.  *2»,X,8„*\.8;, 


it’ 


(») 


where  the  coefficients  ai~  aao  involving  £  and  fx  are  listed  in  Appendix  J. 
When  £-1  the  expressions  for  ai~  aaQ  siaplify  greatly  and  reduce  to  the 
coefficients  Cx-C1#  already  used  in  Ref  [15).  For  r  and  a>3,  it  is  best  to 
carry  out  the  tedious  enuaeration  by  a  computer,  and  a  sample  listing  of  such 
computer  enuaeration  for  £-1.2  and  1  is  presented  in  Appendix  K. 

B.  Clamped  Plate 

For  the  clamped  plate,  using  T*s*t*fc^ (£4+$£*/3-»-l)/3ott^(l-**i )(£*+!)  as  the 
T*.  the  modal  equations  obtained  from  Eqs.  (3.19-3.21)  and  (3.23)  have  the 
dimensionless  fora  in  parallel  to  Bq.  (3) 


3*W 


(CP-1) 


a  +  ~Sr  prs  +  ^i^mn^rs 


(CP-2)  (Cp-3)  (CP-4) 


*  3^(lV)Jr,(»„. 


(CP -5a)  (CP-5b) 

-  OS1+2&LI .9* l ) j  (|f  )+£  (n  \i 

(^a+1j  olp  sl  nr  *s'  mn,Ji 

(CP-6) 


‘rs 


-  ll2i}£lll+yZ&:ll  y  {w  Tpq  ) 

(£*+l )  r#l  “n'  (£«pW) 

(CP-7) 
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lela&siiim  (e  )j  , 

“V'W'rs 

(CP-8) 


(10) 


Because  of  the  parallelism,  Eq.  (10)  has  exactly  eight  terms  CP  (clamped 
plate)  1-8  with  the  same  physical  interpretation  as  SSP  1-8. 

In  reality,  however,  the  enumeration  of  Eq.  (10)  would  generate  many  more 
terms  than  Eq.  (3)  because  the  terms  CP-4,  CP-6,  and  CP-8  are  non-diagonal. 

In  any  event,  the  lowest-order  components  for  r  and  s-1  or  3  are  given  by 

^rs-  Pr»*  »,W  (CP“S)r8“  Vra 


(^♦l) 

_  HzulgL ieh 22I/3+11  9  (_Iesl_)  .  i&zjL/z+iirr,  1  .  0 

(^+1)  r8VpW  StfM)  1  aJrs  ‘ 


(11) 


First,  the  components  of  CP-4  are 


-4s)*,,*  -r>»: 


33 


«A,"  a#1***  f1"..*  sV.dw,,- 

fi.fV*  §*>«„ 


(12) 


Second,  the  components  of  CP-6  are 


Id4*, ^»sl*  fUd4*!)*,,- 


as 


(13) 


Third,  the  components  of  CP-8  are  also  given  by  Eq.  (13)  but  with  replaced 


by  8.,,  arc 
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Fourth,  the  coaponents  CP-7  are 

*.".1*  2t“.»  • 

*.’V  *."«• 

*A,*  e.o"..  • 

where  g’s  are  given  by  Bq.  (F8)  of  Appendix  F.  For  instance,  we  see  that 

gt —  14.222(T0#+Ta0/»"*)/4  ♦  14.222(T()4+T40/J'a  )/16  +  14.222Taa/4(0*+l) 

-  7 . 1 11 <Ta4 (0*+4 ) ~l +T4a (40*+l)_i}/4, 


(14) 


(15) 


(16) 


and  the  remaining  g's  are  expressed  similarly.  Lastly,  instead  of  the 
analytical  expressions,  we  present  here  the  computer  listing  of  CP-5  with  the 
coefficients  evaluated  for  0-1.2  and  fia*0.1 


( CP-5 )  «  0 . 76804E+02xW( 1 , 1 )V ( 1 , 1 )N ( 1 , 1 ) -0 . 14892E+03xW( 1 , 1 )W( 1 , 1 )W( 1 , 3 ) 

11  -0. 18330E+03xW(l, 1)0(1, 1)#(3, l)+0.43305E+02x*f(l,l)W(l,l)W(3, 3) 
+0 . 38925E+03xtf(  1 , 1 ) W( 1 , 3 ) W( 1 , 3 ) +0 . 33102B+03X# ( 1,1)N(1,3)N(3,1) 
-0 . 28888E+03xW( 1 , 1 )lf  ( 1 , 3 )W( 3 , 3 ) +0 . 55192E+03XW ( 1 , 1 )W( 3 , 1 )W( 3 , 1 ) 
-0 . 36109E+03xW( 1 , 1 )  W(  3 , 1 )  W(  3 , 3 )  +0 . 46008E+03xW( 1 , 1 )W( 3 , 3 ) W( 3 , 3 ) 

-0 . 17082E+03XW ( 1 , 3 )W( 1 , 3 )H ( 1 , 3) -0 . 42937E+03XH ( 1 , 3 )W ( 1 , 3 ) W( 3 , 1 ) 

+0 . 17620E+03xW( 1 , 3 )  W(  1 , 3  )W  (  3 . 3 )  -0 . 45437E+03XW ( 1 . 3 )K ( 3 , 1 ) W( 3 , 1 ) 
+0 . 88712E+03xW( 1 , 3 )W( 3 , 1 )W( 3 , 3 ) -0 . 27129E+03 XW( 1 , 3 )  W(  3 , 3  )tf  ( 3 . 3 ) 
-0 . 32645E  +03xH  ( 3 , 1 )  W(  3 , 1  )W  ( 3 , 1 )  +0 . 35047B+03xW(  3 , 1 )W( 3 , 1 )M( 3 . 3 ) 
-0.54223B+03xM(3,l)W(3.3)W(3,3)-0.92581E+01xW(3,3)lf(3,3)W(3.3), 


(CP-5)1#  — 0 .49640B+02xW( 1 , 1 )¥( 1 , 1 )W( 1 , 1 ) +0 . 38925E+03xW( 1 , 1 )W( 1 , 1 )W( 1 , 8) 
1#  +0 . 16551E+03XW  ( 1 . 1  )W(  1 , 1  )W(  3 , 1 )  -0 . 14444E+03x*f(  1 . 1  )W  ( 1 , 1  )W(3 , 3) 
-0 . 51247B+03xtf(  1 , 1  )W(  1 , 3  )W(  1 , 3 )  -0 . 85874E+03xW(  1 , 1  )M(  1 , 3  )N{  3 , 1 ) 
+0 . 35241E+03xW( 1 , 1  )W( 1 . 3 )W( 3 , 3 ) -0 . 45437E+03X* ( 1 , 1 )*( 3 . 1 )tf  <  3 , 1 ) 
+0.88712B+03xW(l,l)W(3,l)W(3,3)-0.27129B+03xlf(l,l)W(3,3)M<3.3) 
+0 . 14310E+04xW( 1 , 3 )W( 1 , 3 )* ( 1 , 3 ) +0 . 33868E+03X*( 1 , 3 )W ( 1 , 3 )W( 3 , 1 ) 
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-  0 . 16894E+04xW( 1 , 3 )W( 1 , 3 ) W( 3 , 3 ) +0 . 17030E+04xW( 1 , 3 )W(3 , 1 )t»(3 , 1 ) 

-0 .78637E+03xW( 1 , 3)W{3, 1 )W(3 , 3) +0. 31710E+04xW( 1 ,3)W(3 ,3)W(3 , 3) 

+0 . 18257E+03xW(3 ,  1 )  W(  3 , 1 )  W(  3 , 1 )  -0 . 75391E+03xW( 3 , 1 )W( 3 , 1 )W ( 3 , 3 ) 

-0 . 28530E+02xW( 3 ,  1 )  W(  3 , 3 )  W(  3 , 3 )  -0 . 3231 1E+03xW( 3 , 3 )W( 3 , 3 )W( 3 , 3 ) , 

(CP-5)  ^=-0 . 61 101E+02xW( 1,1)W(1,1)W(1,1 )+° . 16551E+03xW( 1,1)W(1,1)W(1,3) 

+0 . 55192E+03xW( 1 , 1 )W ( 1 , 1 )W( 3 . 1 )-0 . 18055E+03xW( 1 . 1  )W(  1 . 1  )tf  (  3 . 3 ) 

-0 . 42937E+03xW( 1 , 1 ) W( 1 , 3 )W( 1 . 3 ) -0 . 90874E+03xW( 1 , 1 )*f  ( 1 , 3) *( 3 , 1 ) 

+0 . 88712E+03xW(  1 , 1 )  W(  1 , 3 )  W(  3 , 3 )  -0 . 97935E+03xW(  1 , 1  )*f(  3 , 1 )  W{  3 , 1 ) 

+0 . 70094E+03xW( 1 , 1  )W(  3 . 1  )W(  3 , 3)  -0 . 54223E+03xW( 1 , 1 ) W( 3 , 3) W( 3 , 3) 

+0 . 11289E+03xW( 1 . 3 )tf  ( 1 , 3 )W( 1 , 3 )+0 . 17030E+04xW( 1 , 3 ) W( 1 , 3 ) W( 3 , 1 ) 
-0.39318E+03xW(  1 , 3)W( 1 ,3)W{3 , 3)+0. 54770E+03xW( 1 ,3)W(3, 1 )W(3 , 1 ) 

-0 . 15078E+04xW( 1 , 3  )W(  3 , 1  )W(  3 , 3 )  -0 . 28530E+02xW( 1 , 3 )W( 3 , 3 ) W( 3 , 3 ) 

+0 . 27674E+04xW(3 .  1  )*f  ( 3 , 1 )  W(  3 , 1 )  -0 . 30I67B+04xlf  ( 3 , 1 )W( 3 , 1 )W( 3 , 3 ) 

+0 . 52367E+04xW( 3 , 1  )K(  3 , 3 )  W(  3 , 3 )  -0 . 13402E+03X* ( 3 , 3 )W( 3 , 3 ) W( 3 . 3 ) , 

( CP-5 )#> «  0 . 14435E+02xW( 1 . 1  )*f  ( 1 , 1 )  W(  1 , 1 )  -0 . 14444B+03xW( 1 , 1 )W ( 1 , 1 )*f ( 1 , 3 ) 

-0 . 18055E+03xW( 1 , 1 ) W( 1 , 1 )W( 3 , 1 ) +0 . 46008E+03xW( 1 , 1 )W ( 1 , 1 )W(3 .3) 

+0 . 17620E+03xW( 1 , 1 )W( 1 ,3)W ( 1 , 3) +0 . 88712B+03xW( 1 , 1  )W( 1 . 3)W( 3 , 1 ) 

-0 . 54259E+03xW( 1 , 1 ) W( 1 , 3 ) W( 3 , 3 ) +0 . 35047E+03xW( 1 , 1 ) W( 3 , 1 ) W( 3 , 1 ) 

-0 . 10845E+04xW( 1 . 1  )W(3 , 1  )*f  ( 3 , 3)  -0 . 27774E+02xW( 1 , 1 )*f ( 3 , 3)W(3 , 3) 

-0 . 56313E+03xM( 1 , 3 ) W( 1 , 3 ) W( 1 , 3 ) -0 . 39318E+03xW( 1 , 3 )W ( 1 , 3 )W( 3 , 1 ) 

+0 . 317l0E+04xW( 1 , 3 )  W(  1 , 3 )  W(  3 , 3 )  -0 . 75391B+03xW( 1 , 3 )W ( 3 , 1 ) W( 3 , 1 ) 
-0.57060E+02xW(l ,3)W(3, 1 )W(3,3)-0.96934E+03xW(l ,3)W(3,3)Vf(3,3) 

-0 . 10062E+04xW( 3 . 1 ) W( 3 . 1 ) W(3 . 1 ) +0 . 52367E+04xW(  3 , 1 ) W( 3 . 1 )W( 3 , 3 ) 
-0.40205E+03xW(3,l)W(3,3)W(3.3)+0.42215E+04xW(3,3)lt(3,3)W(3f3).  (17) 

This  therefore  completes  specification  of  Eq.  (11)  for  r  and  s«l  or  3,  which 

clearly  involves  a  lot  aore  terms  than  Eq.  (4)  for  the  simply-supported  plate. 
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V.  Hamiltonians  for  The  Kodal  Equations 

Even  for  only  four  components  of  Wr8,  the  modal  equations  have  a 
bewildering  number  of  constant,  linear,  and  cubic  terms.  Hence,  one  cannot 
help  but  wonder  if  the  modal  equations  have  included  only  those  terns  which 
rightfully  belong  to  them  and  left  out  none  inadvertently.  The  purpose  of 
this  section  is  to  show  the  internal  consistency  of  modal  equations.  We  do 
this  by  way  of  constructing  the  Hamiltonians  for  simply-supported  and  clamped 
plates,  from  which  Eqs.  (4.4)  and  (4.11)  can  be  recovered  by  the  Hamiltonian 
equations  of  motion.  Therefore,  the  readers  who  have  sufficient  faith  in  the 
accuracy  of  Sec.  IV  may  skip  this  section  without  loss  in  continuity. 

A.  Simply-Supported  Plate 

We  define  the  column  vectors  II«(W  •*«,'*»,  .W„),  >®,, )  • 

Mtrice. 


where  a^j  and  b^j  are  given  by  Eqs.  (4.5)  and  (4.6),  and  a  symmetric  matrix 


where  hj  are  given  by  Eq.  (4.8).  With  the  4x4  unit  matrix  I  Eq.  (4.4)  can 
readily  be  put  in  matrix  form 


W  +  7fl«  -  IF  +  S4W  ♦  (SSP-5)rg  -  (f*+l)T0S#» 
-  -  »• 


except  for  the  term  SSP-5.  In  fact,  one  can  also  express  SSP-5  in  matrix 


<SSP-5)r3.3*r*{2S„-<l-<.*)S7t— _f 


(^ap*+q8/Pa)a 


Here,  S8a  is  a  4x4  diagonal  matrix  with  the  elements 

{[(^♦2M+*‘a)Wj1+(^+lOp+90~a)Wj#+(908+lOM-*~a)Wj1+9(*a+2|U+0‘a)Wj8I, 
[(*a+10M.rI#'a  )Wj1  +  (^a+l%i+8V“a  )Wj#  +  (9?a+8^+9p'a  )f$t+9(0a+iqu+8p'*  )Wj#  ]  , 
I  ( )W*  +(9fi*+82u+9fi~* )wf  +(810a+18 )Wf  +9(90a+l<Wa )Wf.  )  , 
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9[  (0*+2 )*^  j +{^a+10/i+9/J"a  )Wj#  +  ( Sf?+10n+fi~z  +9(/Ja+2ji+/?'a  )Wj#  ]}- . 


And,  S7[.J  denotes  a  symmetric  Matrix  with  the  element  h^. 

In  the  absence  of  viscous  damping,  it  is  straightforward  to  deduce  the 
Haailtonian  for  Eq.  (1).  To  do  this  requires  defining  a  column  vector  Q  with 
the  components  qt*Wia,  q^-W^,  q^-W^,  q4“WJ3,  and  a  column  vector  P  with  the 
conjugate  variables  P1»fc11,  pa-Wi3,  P3-W3Jl,  P4**,*»-  Ignoring  the  SSP-2,  the 
Haailtonian  H_  of  the  simply-supported  plate  is  given  by 

s 


v  KpTlP  +  fiTs*Q  "  (^+1  )t  ots  o  -  U-.tfJb*(g?.+i}  qts  (  Tpq  )B\ 
8  451  4  0  8  4  7  0*pa+  q8  * 


~  fiTlF  -  ^r1  “V  +  V 


(3) 


Here,  hg  denotes  contribution  by  the  SSP-5,  to  be  determined  presently. 
Although  it  is  tempting  to  formally  define  h_  by  preaultiplying  Eq.  (2)  by  HT, 
the  resulting  quartic  expression,  though  correct  in  form,  does  not  have  the 
correct  coefficients.  This  is  because  the  quartic  form,  unlike  a  quadratic 
fora,  generate  nonunifora  numerical  factors  under  differentiation.  We  shall 
therefore  deduce  h„  directly  from  Eq.  (4.9)  [15] 

ha  ♦“»«  ♦“,«  "AV. 

'\A%%  sw*,  *\A  *\A< 

*",.^>'4  **,.44  **,,151.  *».!,<•  <4> 

It  is  siaple  to  checked  that  Eq.  (4.4)  is  recovered  from  the  Haailtonian 
equations  of  motion 

qj-SH/dPj,  Pj=  -3H/3qj,  (5) 


when  H»H8  is  given  by  Eqs.  (3-4). 

B.  Clamped  Plate 

Similarly,  Eq.  (4.11)  can  also  be  put  in  matrix  fora 


W  +  '*IW  -  IF  4  C.N  +  (CP-5)  - 

4  r*  (*M) 

_  (  Ie3— )H  _  Il-t%?f/3+li 

(0*+i)  7  p j^+q*  3(p*+i) 
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where  gj  are  given  by  Eq.  (4.15).  Hence,  the  Hasiltonlan  Hc  for  the  clasped 
plate  becomes 

H  4/pTIP  +  0TC  fi  -  ll£!±2fifZ3illT  0TC  D  .  n-elguhz&atnatc  (  -Tpg — )q\ 

c^lrxr  "V*  {p*+1)  VV*  (^1)  S  Vl*+«r  ' 


fiTIF  -  ♦  h£ 

3(0*+l)  *  c 


(7) 


And,  the  contribution  by  CP-S  deduced  frow  Eq.  (4.17)  is 


hc«  0. 19201E+02xq(l)q(l)q(l)q(l) 
-0.61101E+02xq(l)q(l)q'l)q(3) 
+0 . 19463E+03xq( 1 )q( 1 )q(2)q( 2) 
-0. 14444E+03xq(? )q(l )q(2)q(4) 
-0.18055E+03xq(l)q(l)q(3)q(4) 
-0.17082E+03xq(l)q(2)q(2)q(2) 
+0.17620E+03xq(l)q(2)q(2)q(4) 
+0.88712E+03xq(l)q(2)q(3)q(4) 
-0.32645E+03xq(l)q(3)q(3)q(3) 
-0.54223B+03xq(l)q(3)q(4)q(4) 
+0.35776E+03xq(2)q(2)q(2)q(2) 
-0.56313E+03xq(2)q(2)q(2)q(4) 
-0.39318F+03xq(2)q(2)q(3)q(4) 
+0.18257B+03xq(2)q(3)q(3)q(3) 


-0 . 49640E+02xq(l )q( 1 )q( 1 )q(2) 
+0.14435E+02xq(l)q(l)q(l)q(4) 
+0. l6351E+03xq( l)q(l)q(2)q(3) 
+0. 27596E+05xq(l )q{ 1 )q(3)q(3) 
+0.23004E+03xq(l)q(l)q(4)q(4) 
-G.42937E+03xq( 1 )q(2)q(2)q(3) 
-0.45437E+03xq( l)q(2)q(3)q(3) 
-0 . 27129E+03xq( 1 )q( 2 )q( 4 )q(4 ) 
+0.35047B+03xq(l)q(3)q(3)q(4) 
-0.92581E+01xq( 1 )q(4)q(4)q(4) 
+0. 11289E+03xq(2)q(2)q(2)q(3) 
+0.85149E+03xq(2)q(2)q(3)q(3) 
+0. 1585SE+04xq(2)q(2)q(4)q(4) 
-0.75391E+03xq(2)q(3)q(3)q(4) 
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-0.28530E+02xq(2)q(3)q(4)q(4)  -0.32311E+03xq(2)q(4)q(4)q(4) 
+0.69186E+03xq(3)q(3)q(3)q(3)  -0.10062B+04xq(3)q(3)q(3)q(4) 

+0. 26183E+C4xq(3)q(3)q(4)q(4)  -0.13402E+03xq(3)q(4)q(4)q(4) 

+0. 10554E+G4xq(4)q(4)q(4)q(4) .  (3} 

Again,  the  Haailtonian  equations  of  aotion  rederive  Eq.  (4.11)  when  H*H„. 

c 

Suaaing  up,  the  existence  of  H_  and  H_  strongly  suggests  the  internal 
consistency  of  aodal  equations  for  the  siaply-supported  and  claaped  plates, 
whereby  no  teras  have  been  added  or  left  out  inadvertently.  Unfortunately, 
this  does  not  prove  the  absolute  correctness  of  the  aodal  equations  derived  in 
Sec.  III.  The  reason  is  that  the  Haailtonians  are  deduced,  in  part,  froa  the 
aodal  equations.  Hence,  errors  that  are  consistent  with  the  Haailtonian 
foraulation  can  reaain  undetected. 
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VI.  Mid-Plate  Teaperature  Variation  and  Gradient 

To  proceed  further  with  the  nodal  equations,  it  is  necessary  to  specify 
Tm  and  0r#  which  appear  in  Eqs.  (4.4)  and  (4.11).  Returning  to  Eq.  (2.3), 
we  first  split  T(x,y)  into  the  unifora  teaperature  tQ  and  teaperature  varia¬ 
tion  ty(x,y)  over  the  plate 

T  -  tfl  +  tv(x.y).  (1) 

as  in  Eq.  (3.2).  We  then  convert  $  to  the  corresponding  teaperature  differen¬ 
tial  9  (see,  Eq.  (4.2)) 

0  *  h0.  (2) 

In  view  of  Eqs.  (1-2),  we  nay  pot  Eq.  (2.3)  in  the  following  fora 

T  -  t0  +  ty(x,y)  ♦  ®(x.y),  (3) 

where  Z*z/h  ranges  over  (1/2, -1/2).  Since  tQ  will  be  assuaed  nonzero,  we  use 
it  as  the  aain  thermal  paraaeter  to  express  the  aagnitudes  of  ty  and  ® 

ty- «vt#fv(x,y),  0  -  «gt0fg(x,y).  (4) 

Here,  the  scaling  factor  4y  defines  the  aagnitude  4ytQ  of  teaperature  varia¬ 
tion  whose  profile  over  the  plate  is  given  by  fy(x,y).  Similarly,  4  tQ  is  the 
magnitude  of  temperature  gradient  through  the  plate  thickness  and  fg(x,y) 
represents  its  distribution  over  the  plate.^  Upon  introducing  Eq.  (4)  Into  Eq. 
(3),  we  have  the  diaensiunless  teaperature  T  -T/T* 

T  -  T0  +  («vT0)fy(x.y)  +  Z(4gT0)fg(x,y),  (5) 

where  T*  is  T*  for  a  siaply-supported  plate  and  T*  for  a  claaped  plate. 

First  of  all,  for  the  teaperature  variation  over  the  plate,  we  find  by 
coaparing  Eq.  (9)  with  Bq.  (3.2)  that 


(4wTn)f„(x.y)«  ST  if*  T  cospnx  cosqrry, 
0  p-o  q-o  M 


(6) 


which  holds  for  both  the  siaply-supported  and  claaped  plates.  Let  us  examine 
two  examples  of  fy(x,y)  which  are  neither  constant  nor  linear.  As  a  first 
example,  we  let  fy-sinrrxslimy  and  find  from  Eq.  (6)  that 


T 


Si 

pq 


164  T 

_ Jwv*o  _ 

«*(p*-i  >  (*•-!)  (i-«Jm-4j;) 


(p-q-0  excluded) 


(?) 


Now,  for  the  second  example  fy*sinaTrxsin*iry,  the  sum  of  Eq.  (6)  reduces  to  a 
finite  sum  (l/4)(i-cos2sx)(l-cos2iry) .  Hence,  we  have  by  inpaectlon 
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(8) 


T,o-  T0*“  _Taa“  _evTo/4,  Tpq*  0  *for  othcr  p  and  q)* 


Since  the  majority  of  T__  is  zero,  Eq.  (4.8)  reduces  to 

Pm 

x  i* 

(l+p"a),  etc.. 


(9) 


for  the  simply-supported  plate,  and  Eq.  (4.16)  elves  rise  to 

*i-  ^(-t r)t1+*"a+  etc-  <10> 

for  the  claaped  plate  (Note  that  the  coefficient  14.222  in  Eq.  (4.16)  is 
128/9).  And,  the  reaainlng  h's  and  g's  summarized  in  Appendix  L,  all  have 
simplified  expressions. 

Second,  for  the  temperature  gradient  across  the  plate  thickness  we  find 
from  Eqs.  (3.6)  and  (3.18) 

(8_T  )f  (x,y)*  if*  ^enm^tx^ty).  UD 

*  *  »*o  n-o 

where  the  function  nQ  is  i>n  for  a  simply-supported  plate  and  ?n  for  the  clamp¬ 
ed  plate  case.  For  maximal  simplicity,  we  shall  assume  fg*slnrrxsinrry  for  the 
simply-supported  plate  and  obtain 

0tl*  8gTQ/2  an<*  0  (for  other  a  and  n).  (12) 

On  the  other  hand,  fg«sirftrxsii£ny  gives  rise  to 

8.  »  34_T  /8  and  8  *  C  (for  other  m  and  n),  (13) 


for  the  claaped  plate.  Clearly,  8^=0  for  a«n*l  is  generally  not  observed 
under  different  fg.  Using  Eq.  (12)  the  last  term  of  Eq.  (4.4)  has  a  single 
nonzero  component  -(d*+l)a4gTQ/24.  And,  similarly,  Eq.  (4.14)  simplifies  to 


Uta}lt-(**+l)«gT0/2,  Wa ]1# - [Aa J4J if  —4 gT0 //B ,  and  under  Bq.  (13). 

The  relationship  between  the  magnitudes  «yT0  and  4gTfl  is  shown  schema¬ 
tically  in  Fig.  3.  Although  it  appears  at  first  sight  that  the  ty(x,y)  and 
0(x,y)  can  be  assigned  arbitrarily,  this  is  not  the  case  when  they  have  a 
nonunifora  distribution  over  the  plate.  First  of  all,  we  notice  that  a  non- 
uniform  9(x,y)  cannot  exist  unless  ty(x,y)  is  also  nonunifora.  That  is,  4y«0 
implies  6g»0.  Nhat  is  then  the  maximum  value  of  4g,  denoted  by  (4g)aax.  f°r  ® 
given  6y?  This  cannot  be  answered,  in  general,  without  knowing  the  profiles 
fy(x,y)  and  fg(x,y).  However,  in  the  case  of  fy(x,y)«fg(x,y)  (i.e.,  fy«fg» 
8inanxsin®ny  for  the  clamped  plate  discussed),  it  is  readily  seen  that 
(6g)aax^2«v>  Sunakawa  and  Uemura  [13]  have  used  a  parabolic  temperature 
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distribution  for  both  f„  and  f_,  and  consequently  prescribed  (4  „,)„,,■(  4/3  )4„, 
which  iaplles  that  the  upper  surface  teaperature  is  half  the  lower  surface 
teaperature  (see,  Bq.  (42)  in  Ref  [13]). 

In  view  of  the  recent  atteapts  [17,18]  to  generate  various  teaperature 
profiles  by  radiant  heating,  teaperature  variation  profiles  of  the  sort  that 
we  have  considered  here  do  not  appear  at  all  unrealistic.  On  the  other  hand, 
to  the  best  of  our  knowledge  no  atteapt  has  yet  been  aade  to  either  aeasure  or 
iapose  certain  teaperature  gradient  distribution  through  the  plate  thickness. 
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VII.  Single-Mode  Equation  for 

Let  us  examine  the  simplest  case  of  modal  equation  for  Wtl  when  all  other 
Wr8  are  absent.  With  Eqs.  (6.9)  and  (6.12),  Eq.  (4.4)  for  the  simply- 
supported  plate  becomes  after  letting  r-s*l 

where  at  (/*.£)*  |((l-Ma)(0*+l)  +  2{p*+l+2np*)]  as  shown  in  Appendix  J. 
Similarly,  with  the  use  of  Eqs.  (6.10)  and  (6.13),  Eq.  (4.11)  for  the  clamped 
plate  reduces  to 

l}w„ 

+  J(^+2^»/3-H)4gT0]-  0,  (2) 

where 

1  4  l  9L®  [0+g  *)*  (p+4P  4)a  (4 fi+fi  4)*JJ 

Here,  (128/9)d1  (✓0.1,1.2)*76.804  is  the  numerical  coefficient  for  W(l,i;*-tern 
in  (CP-5)lt  as  given  by  Eq.  (4.17).  The  eight  terms  in  Eqs.  (4.4)  and  (4.11) 
have  been  regrouped  into  five  in  Eqs.  (1-2).  In  particular,  note  that  the 
third  term  of  Eq.  (1)  is  the  combined  stiffness  which  subsumes  structural 
stiffness  (SSP-4) ,  global  thermal  expansion  by  uniform  temperature  (SSP-6), 
and  local  thermal  expansion  by  temperature  variation  (SSP-7).  And,  the  last 
term  in  Eq.  (1)  is  the  combined  forcing  of  both  the  applied  external  pressure 
(SSP-3)  and  thermal  moment  owing  to  temperature  gradient  across  the  plate 
(SSP-8).  That  the  temperature  gradient  plays  the  role  of  external  forcing  has 
already  been  observed  by  Boley  and  Weiner  [7]  and  Sunakawa  and  Ueaura  [13]. 

In  a  parallel  fashion,  the  combined  stiffness  and  forcing  terns  in  Eq.  (2)  are 
giver,  by  the  regrouping  of  terms  in  Eq.  (4.11). 

Since  Eqs.  (1-2)  are  qualitatively  similar,  we  may  present  them  in  a 
prototype  form  by  denoting  q«W  and  f»P 

*q  +  «0*q  +  «J(l-s)q  +  wq*  »  f0  +  f,  (3) 

where  «o«(0*+l),  s«T0 [1+  -g(l-*»)«v] .  *-4a1,  f0«  for  the  simply- 

supported  plate;  and  «*« -i|(/»4+8P*/3+l),  S«T0[1»  •g(l-y)6v[l+  ]. 

*»  T^d4*  fQ*  ^(^4+2)»*/3+l)6gT0  for  the  clamped  plate.  Note  further  that  r-«0 
has  been  introduced.  The  combined  stiffness  «*(!-*)  remains  positive  as  long 
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as  the  theraal  loading  is  weak  (pre-buckling  for  s<l),  whereas  it  becoaas 
negative  under  a  strong  theraal  loading  (poet-buckling  for  s>l). 

The  Haailtonian  for  Bq.  (3)  including  the  kinetic  energy  and  potential 
(strain)  energy  can  be  written  down  froa  Eqs.  (5.3)  and  (5.7) 

H  •  (4) 

As  an  illustration,  we  have  shown  in  Fig.  4a  the  potential  energy 
«4-:  (l-s)<f +^q4  of  a  siaply-supported  plate  for  0*1  andp?*0.1,  while  incre- 
aenting  s  froa  0  to  3.  More  visually,  however,  Fig.  4b  depicts  the  potential 
energy  surface  which  is  a  single  well  potential  for  s<l,  but  develops 
syaaetric  double-well  potential  as  s  beccaes  larger  than  unity.  Both  T0  and 
4yTQ  contribute  to  the  paraaeter  s.  We  have  shown  in  Fig.  5  the  threshold 
boundary  of  theraal  buckling  (s*l),  which  intersects  the  T0~axls  at  TQ*1  and 
approaches  the  4y-axis  asyaptotically. 

A.  Theraal ly  Buckled  Nodal  Aaplitude 

For  the  static  problea,  we  retain  only  the  teras  which  do  not  Involve  tlae 
differentiation 

*0*  +  *{(l-s)Q  -  g  -  0,  (5) 

where  g»f0+f  is  the  coabined  forcing.  Under  a  weak  theraal  loading,  the 
coabined  stiffness  is  positive,  hence  only  one  root  of  Eq.  (5) 


is  real.  Note  that  g  aust  be  nonzero  in  Bq.  (6);  otherwise.  Qt*0  is  the  only 
real  root.  Because  of  g*fQ+f,  the  teaperature  gradient  alone  is  sufficient  to 
sustain  theraal  buckling,  even  when  there  is  no  external  pressure  Japosed.  In 
contrast,  for  a  strong  theraal  loading  the  coabined  stiffness  becoaes  nega¬ 
tive,  hence  one  aay  either  include  [16]  or  exclude  [6]  the  coabined  forcing  In 
coaputing  the  buckled  aaplitude.  In  the  latter  case,  the  theraally  buckled 
aodal  aaplitude  is 

V  (7) 


B.  Oynaaic  Considerations 

In  the  pre-buckled  state,  qc«0  is  the  equlllbrlua  state  of  Eq.  (3),  and 
corresponds  to  the  location  of  single-well  potential  energy  in  Fig.  4.  Note 
u»at  Eq.  (3)  is  a  Duffing  equation  (19).  On  the  other  hand,  q0«0  becoaes 
unstable  in  the  post-buckled  state.  Hence,  the  stable  equlllbrlua  states  are 
now  given  by  q^O^ ,  which  corresponas  to  the  locations  of  double-well 
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potential  energy  (Fig.  4).  We  then  rewrite  Eq.  (3) 


*q  +  -  «*(s-i)q  +  *q*  -  f0+  <8) 

which  is  the  so-called  buckle d  been  equation  originally  investigated  by  Holmes 
[20,21],  the  trajectory  of  which  wanders  in  and  out  of  the  potential  energy 
wells  in  a  chaotic  fashion.  By  the  change  of  variable 


we  obtain  from  Eq.  (8) 


dm 

q  -  + 


(9) 


q  +  «0$q  +  2«*(s-l)q  ±  ^(^q*  +  wq*  -  fQ+  f,  (10) 

which  now  represents  oscillation  about  Q g .  Coaparing  the  linear  stiffness 
terms  of  Eqs.  (3)  and  (10),  Schneider  [3]  has  concluded  that  the  natural 
frequency  increases  by  the  factor  /Z  after  a  thermal  buckling.  The  dynamical 
behavior  of  Eqs.  (3)  and  (8)  with  respect  to  s  will  be  investigated  elsewhere. 
At  present,  however,  the  wore  compelling  need  is  to  estimate  the  stochastic 
response,  for  acoustic  excitations  will  be  used  in  the  high-temperature  sonic 
fatigue  test  facility  being  constructed  at  the  Structual  Dynamics  Branch 
(WL/FIBG) . 
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VIII.  Response  Estimation  by  the  Equivalent  Linearization  Technique 

Since  fQ  is  assuned  constant,  the  combined  forcing  is  a  Gaussian  process 
superposed  on  a  nonzero-nean  level  when  f(t)  is  aean-zero  Gaussian.  The 
equivalent  linearization  as  originally  foraulated  by  Caughey  [22]  and  Booten 
[23]  was  built  on  the  zeri-aean  Gaussian  excitation,  Hence,  it  cannot  be 
applied  directly  to  Eqs.  (7.3)  and  (7.10)  when  fQ*0.  In  Sec.  Villa,  we  re¬ 
strict  ourselves  to  fo»0  and  apply  the  equivalent  linearization  in  its  faail- 
iar  fora.  The  case  of  fQ?£0  will  be  presented  in  Sec.  VUIb,  after  extending 
the  equivalent  linearization  technique  to  Gaussian  forcing  with  a  nonzero  aean 
(see.  Appendix  N). 


A.  Acoustic  Loading 

For  s<l  let  us  denote  by  qjjn  the  aaplltude  of  linearized  Eq.  (7.3). 
Under  fQ»0  the  aean  square  aaplltude  is  given  by  Eq.  (N12)  of  Appendix  N 


<"5m> 


(i) 


where  <  >  is  the  statistical  (tine)  average.  Here,  t^(u)  is  the  power 
spectral  density  of  Gaussian  forcing  expressed  in  angular  frequency  «  [8] . 
Assuaing  that  is  more  or  less  flat  around  the  sharp  resonance  peaks  at 

u*±uQJi-s ,  Eq.  (1)  has  the  alternate  fora 


% 


gff 

2*«J(l-s) 


(2) 


using  (w)"Eff(f )/2rt ,  where  f  is  frequency.  The  arguaent  of  gff(f)  is 
suppressed,  for  gff  is  assuaed  a  constant  power  spectral  density.  Now,  for 
the  full  nonlinear  Eq.  (7.3)  the  aean  square  aaplltude  can  be  estimated  by  the 
equivalent  linearization  technique 


<qa>* 


1^<0^An> 

oj(l-s) 


(3) 


as  given  by  Eq.  (M18)  of  Appendix  N.  As  expected,  for  a  saall  «,  i.e.,  weak 
nonlinearity,  Eq.  (3)  reduces  to 


«f> » «s5ln>. 

whereas,  when  the  cubic  nonlinearity  is  doainant  we  find  that 


(4) 


<q*>  *  v^«Iiin>.  (5) 

Hence,  <q*>  increases  aore  gradually  as  the  square  root  of  the  aagnitude  of 
forcing  power  spectral  density  in  Eq.  (5),  rather  than  linearly  in  Eq.  (4). 
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As  the  r.ontheraal  (s*0)  reference  case.  Me  have  evaluated  <q*jin>  and  <q*> 
froa  Bqs.  (2-3)  for  £-1,  m*“0.1,  and  f-0.04,  and  presented  in  Fig.  6  the 
aaxiaua  mean  square  displaceaents  coaputed  by 


<(waax/h)8>  3  b<<*lin>  or  b<<5a>' 


(6) 


where  b»4  for  a  siaply-supported  plate  and  b-64/9  for  a  claaped  plate.  In 
the  figure,  the  straight  lines  originating  froa  the  origin  are  the  linear 
input-output  relation  given  by  Eq.  (2).  In  contrast,  the  aean  square  displace- 
■ent  predicted  by  Eq.  (3)  increases  aore  gradually  due  to  the  nonlinear  energy 
sharing  provided  by  the  equivalent  linearization.  For  a  given  g^f,  the  siaply- 
supported  plate  has  a  larger  nean  square  displaceaent  than  the  corresponding 
claaped  plate,  as  already  pointed  out  by  Nei  [11]. 

In  the  cheraal  case  (s>0),  however,  the  aaxiaua  aean  square  displaceaent 
does  Increase  with  s,  as  indicated  by  the  three  values  of  s*0,  0.5,  and  0.9  in 
Fig.  7.  However,  we  notice  that  the  increase  is  aore  pronounced  in  the  claap¬ 
ed  plate  than  in  siaply-supported  plate.  In  both  plate  cases,  the  net  theraal 
contribution  diainshes  as  g^f  becoaes  large,  which  is  supported  by 

<<*>  *  ’  (7) 


obtained  froa  Eqs.  (2-3)  in  the  liait  as  gff  •*».  Eq.  (7)  states  that  the  res¬ 
ponse  is  essentially  Independent  of  theraal  effects  when  acoustic  excitations 
becoae  very  intense. 

Next,  let  us  exaaine  the  pair  of  Eqs.  (7.9-7.10)  for  s>l.  Now,  denoting  by 
qlin  *be  aRPlitude  of  linearized  Eq.  (7.10),  we  have  in  parallel  to  Eq.  (2) 


as 


gff 

4fuJ(s-l) 


(8) 


The  difference  in  nui  cal  factors  of  Eqs.  (2)  and  (8)  is  due  to  the 
factor  associated  with  the  post-buckled  natural  frequency  in  Eq.  (7.10).  It 
is  siaple  to  check  that  the  equivalent  linearization  goes  through  just  as  in 
the  pre-buckled  case,  for  the  square  tera  ±3*<^q*  in  Eq.  (7.10)  has  no  effect 
under  the  Gaussian  assuaption.  Hence,  in  analogy  to  Eq.  (3)  the  aean  square 
aaplitude  is 


Because  of  Eq.  (7.9),  the  total  aean  square  aaplitude  is  sua  of  the  buckled 
plate  position  and  aaplitude  fluctuations 

<q*'  *  <£  ♦  <q*>,  (10) 
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since  Qa  and  q  are  statistically  independent. 

The  total  mean  square  displacement  shown  in  Fig.  8  increases  greatly  with 
s.  Moreover,  the  increase  is  considerably  more  pronouned  in  a  clamped  plate 
(Fig.  8(b))  than  simply-supported  plate  (Fig.  8(a)).  However,  this  rapid  in¬ 
crease  is  mainly  due  to  the  contribution  of  in  Eq.  (10).  Hence,  to  separate 
the  contribution  of  from  that  of  <q*>.  we  presented  in  Fig.  9  the  total 
mean  square  displacement  over  the  range  of  s*(0,  3)  for  a  particular  value  of 
gjf=l.  We  see  from  the  figure  that  not  only  does  (£  increase  linearly  with 
(s-1),  but  also  <qa>  actually  decreases  steadily  in  the  range  of  s>l  from  its 
maximum  value  at  8*1.  The  latter  is  predicted  by  the  asymptotic  form  of  Eq. 
(8),  <q*>  a  1/s,  as  s  ■*».  Hence,  the  mean  square  displacement  falls  off  as 
1/s  in  the  post-buckled  state.  This  has  been  borne  out  by  the  Monte-Carlo 
simulation  of  Choi  and  Vaicaitis  [12],  in  which  the  amplitude  of  displacement 
fluctuations  is  considerably  small  whenever  the  displacement  and  stress  time- 
histories  jump  off  to  thermally  buckled  positions. 

The  composite  view  of  total  mean  square  displacements  is  shown  in  Fig  10 
over  the  forcing  range  of  Kff3(0,  5).  Notice  that  the  clamped  plate  builds  up 
the  total  mean  square  displacement  much  more  rapidly  than  a  simply-supported 
plate.  However,  this  rapid  buildup  in  the  post-buckled  state  is  mostly  due  to 
the  contribution  of  buckled  plate  displacement,  as  depicted  in  Fig.  11. 

B.  Combined  Acoustic  and  Thermal  Loading 

The  analysis  of  Sec.  Villa  was  restricted  to  f0*0.  Clearly,  this  is  valid 
when  a  plate  is  maintained  in  uniform  temperature  or  there  is  no  heat  flux 
through  the  plate  thickness.  The  latter  is  possible  when  one  side  of  the  plate 
is  heated  while  the  other  side  is  insulated.  However,  under  a  nonuniform  heat 
flux,  the  temperature  gradient  will  have  a  certain  distribution  over  the 
plate,  which  is  neither  constant  nor  linear.  To  incorporate  fo  into  the  acous¬ 
tic  loading,  one  must  extend  the  equivalent  linearization  to  the  nonzero-mean 
Gaussian  excitations.  Let  us  denote  by  qQ  the  amplitude  under  the  combined 
forcing  fQ+f(t)  and  split  it  into  two  parts 

qc  *  x  +  y.  (11) 

The  implicit  assumption  is  that  x  is  the  response  to  fQ  and  y  to  f(t).  We 

first  replace  q  in  Eq.  (7.3)  by  q  and  then  Introduce  (11)  to  obtain 

c 

“y  +  w0fy  +  {wo(l-*)+®c'^ }y  +  Sexy8  +  jey3  +  w2(l-s)x  ♦  k1?  -  f  «  f,  (12) 
Applying  the  equivalent  linearization  technique  to  Eq.  (12)  gives  rise  to 
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(13) 


<y*> 


12K<<](lln> 

^(1-3) 


given  by  Bq.  (N17)  of  Appendix  N.  Furthermore,  the  x  In  Eq.  (13)  is  given  by 
Eq.  (N20) 


The  pair  of  Eqs.  (13)  and  (14)  can  be  solved  simultaneously  for  <f>  and  x. 
Afterwards,  the  total  mean  square  of  qQ  is  computed  by  sum  of  the  squared  x 
and  mean  square  fluctuations  owing  to  acoustic  loading 


«£>  *  x2  +  <y*>.  (15) 

This  is  the  amplitude  response  to  the  combined  acoustic-thermal  forcing.  Note 
that  when  fQ-0  Eq.  (14)  becomes  x-0.  Then,  Bqs.  (13)  and  (15)  degenerate  to 
Eq.  (7.13)  which  was  derived  in  Sec.  Villa  under  the  assumption  of  fQ*0. 

Recall  that  in  Sec.  VII  both  TQ  and  6yT0  are  represented  collectively  by  a 
single  parameter  s.  Under  the  combined  forcing,  however,  one  must  also  specify 
the  magnitude  of  temperature  gradient  6gTQ ,  in  addition  to  s.  For  simplicity, 
we  shall  assume  4e*<5v  in  what  follows.  For  three  values  of  TQ=  0.2,  0.5,  and 
0.8,  we  choose  <5y  large  enough  for  s  to  be  about  0.96  and  compute  the  total 
mean  square  amplitude  by  Eqs.  (13-15).  The  results  are  summarized  in  Pigs.  12 
and  13  for  the  simply-supported  and  clamped  plates,  respectively.  Three  plots 
are  presented  in  each  frame  of  the  figures.  The  first  two  plots,  referring  to 
the  left  ordinate,  are  the  maximum  mean  square  displacement  b<q^>  under  tne 
combined  loading  (solid  curve)  and  b<q*>  under  the  acoustic  loading  alone 
(dotted  curve).  Here.  b*4  for  the  simply- supported  plate  and  b»64/9  for  the 
clamped  plate.  One  observes  that  the  thermal  loading  is  generally  unimportant, 
unless  acoustic  excitations  are  very  weak.  Now,  to  accentuate  the  difference 
of  the  first  two  plots,  the  third  plot  with  respect  to  the  right  ordinate, 
shows  the  difference  of  the  first  two  in  percent,  i.e., 

A<(Wmax/h)a>*100(<q^>-<qa>)/<q|>.  All  the  third  plots  fall  off  rapidly  from 
the  initial  100k  as  gff  increases.  Even  in  the  most  thermal-loading  sensitive 
case  (Fig.  13a),  we  find  that  A<(wBax/h)*>  drops  down  to  4*  when  g^f 
approaches  unity. 

C.  RMS  stress  and  strain 

The  lesson  of  Sec  VUIb  is  that  the  temperature  gradient  across  the  plate 
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which  appears  as  an  additional  loading,  is  generally  unimportant,  as  long  as 
the  acoustic  loading  is  dominant.  Quantitatively,  for  gff>l  me  can  ignore 
the  effect  of  temperature  gradient,  and  hence  consider  only  the  uniform  tem¬ 
perature  and  temperature  variation  over  the  plate.  Nevertheless,  since  the  «T 
tern  appears  explicitly  in  Eqs.  (2. 1-2. 2),  one  suspects  that  the  temperature 
gradient  night  contribute  to  the  stress  and  strain  tensors  as  equally  as  the 
uniform  temperature  and  temperature  variation.  The  purpose  of  this  subsection 
is  to  substantiate  this  suspicion. 

To  this  end,  we  express  the  normal  stress  and  strain  components  involving 
only  the  W4l  in  the  following  form 


°x]  „  t£e£  fexl 

°yj  *  KJ  ’ 

'cxl  ,  «¥  f*x1 

L*yJ  #  KJ’ 


(16) 

(17) 


where  the  dimensionless  cr  ,  o  ,  t  ,  and  t  are  given  by  Eqs.  (012-015)  of 

*  Jr  A  Jf^  ^ 

Appendix  0.  It  should  be  pointed  out  that  ax  and  oy  are  symmetric  in  x  and  y 
(when  one  ignores  the  parameters  p  and  p) ,  and  so  are  *x  and  ey.  We  shall 
therefore  examine  here  only  ox  and  *x  as  a  function  of  x,  under  a  fixed  y-1/2 

Simply-supported  plate 


o_*=- 


(l3-!)  T( 
12(l-*a) 


{l  +(1-M)(^)[l  -  _l_Co92nx] 

+  ■££l^~^(sinnx)q  +  +  p*\if , 

(1-M2)  ^(l-M*)  > 


+ZigSinwx 


(18) 


*x*  i§( T*p )® v^o{ fv“  ^[(It^COsZex) 
+  ^(2 p*+  pcos2nx)qa. 


^COS2rrxl\  +  220*  (silBtx)q 
(**+1)  JJ 

(19) 


Clamped  plate 


Cf 


M 

X 


■f 


ip'+ifl*,  3ri)T  f  f 

ft '  u",(^[1  - 

— ^”-{£*cos27rx  +  psinsnx]q 
3(  !-*»*) 


+  £  +  £  +  POi2nx 
9*-16(l-Ma)  4  2(p*p~l)a 


— | - cos2nxl  +Z<5  »inanx 

U»*+l)  J  B 


CO»2nx 
(p+4p~l )a 


♦ 


-  CQ847TX.A  a 

4<4WM*' 


(20) 
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+  +  i(^*+MC082 nx)  +  ^(P*-MC084rrx)  -  S2S&Z. 

2(P+P  ) 

Eqs.  (18-21)  are  quadratic  expressions  in  q,  which  we  shall  express  in  the 
following  skeleton  form 


IT  =  eo  +  cxq  +  Caq*.  (22) 

where  IT  represents  any  one  of  ax,  oy,  *x,  and  cy  and  CQ ,  Ci ,  and  Ca  are  the 
corresponding  coefficients  involving  x  and  other  parameters.  Under  the  assump¬ 
tion  that  q  is  zero-mean  Gaussian,  the  mean  square  TT  is  given  by 

<n#>  ‘  +  «'f+2C0Ca)<qa>  4  3(^(<q8>)a1  (23) 

in  terms  of  the  <qa>  already  presented  in  Figs.  12  and  13. 

To  be  specific,  let  us  choose  Figs.  12b  and  13b  as  typical  cases  of  the 
simply-supported  and  clamped  plates,  and  also  consider  gff**l»  As  discussed 
in  Sec.  VHIb,  at  this  acoustic  loading  the  effect  of  thermal  loading  is 
negligible  on  the  mean  square  amplitude.  Using  the  parameter  values  of  Figs. 
12b  and  13b,  we  first  compute  the  extreme-fiber  stress  and  strain  by  setting 
Z-l/2  in  Eqs.  (18-21),  then  estimate  the  mean  square  stress  and  strain  by  Eq. 
(23),  and  finally  present  in  Figs.  14  and  IS  the  root-mean-square  stress  and 
strain  distributions  in  the  x  coordinate.  To  quantify  the  role  of  thermal 
terms,  four  plots  are  shown  in  each  frame  of  ’’igs.  14-15-  The  first  (dotted 
curve)  is  the  nonthermal  reference  case  (TQ*  cv*4g*«0);  the  second  (solid  curve 
with  0)  includes  only  the  uniform  temperature  (TQ#0,  4v»4g»0);  the  third 
(solid  curve)  involves  both  the  uniform  temperature  and  temperature  variation 
over  the  plate  (T0*0,  4^0,  4£«0);  and  the  last  (solid  curve  with  X) 
represents  the  fully  thermal  case  (TQ*0,  S^O,  4^*0).  Figs.  14-15  show  that 
three  thermal  terms  do  not  contribute  additlvely  to  the  stress  and  strain 
tensors.  In  other  words,  the  thermal  terms  do  not  necessarily  bring  about 
Increased  stress  and  strain  over  the  entire  x  range.  Moreover,  the  effect  of 
temperature  gradient  which  we  have  written  off  as  unimportant  in  Sec.  VI lib, 
becomes  the  main  contributor  to  the  streso  components. 
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IX.  Assessment  of  Random  Single-Node  Dynamics 

For  a  pre-buckled  plate  (s*l),  the  three  terms  of  uniform  temperature, 
temperature  variation  over  the  plate,  and  temperature  gradient  across  the 
plate,  all  contribute  to  increasing  the  mean  square  displacement,  although  the 
increase  due  to  the  temperature  gradient  is  insignificant  unless  the  acoustic 
loading  is  very  weak.  On  the  other  hand,  the  thermal  terms  do  not  necessarily 
bring  about  increased  stress  and  strain  components,  as  evidenced  by  the  de¬ 
creased  root-mean-square  strain  and  stress  tensors  in  some  region  of  the 
plate.  Moreover,  the  temperature  gradient  contributes  significantly  to  the 
stress  tensor,  while  its  effect  was  negligible  on  the  mean  square  displace¬ 
ment.  The  roles  played  by  three  thermal  terms  are  briefly  summarized  in  Table 
1. 


Table  1.  Summary  of  the  thermal  effects  for  s<l 


Uniform  temperature 

Temperature  variation 
over  the  plate 

Temperature  gradient 
across  the  plate 

Displace¬ 

ment 

In  the  absence  of  temperature  gradient, 
the  uniform  temperature  and  temperature 
variation  are  lum; sd  into  the  parameter 
s  (Fig.  8).  Small  increase  for  a  sim¬ 
ply-supported  plate  and  moderate  in¬ 
crease  for  a  clamped  plate. 

Negligible  for  both 
the  simply-supported 
and  clamped  plates, 
except  men  the  acous¬ 
tic  loading  is  weak 
(Figs.  12  and  13). 

Stress 

Very  small  increase  for  a  simply-sup- 
ported  plate  and  moderate  increase  for 
a  clamped  plate  (Figs.  14  and  15). 

Significant  for  both 
the  simply-supported 
and  clamped  plates. 

Strain 

Moderate  increase 
for  both  the  sim¬ 
ply-supported  and 
clamped  plates. 

Significant  in¬ 
crease  for  both  the 
simply-supported  and 
clamped  plates. 

Negligible  for  both 
the  simply-supported 
and  clamped  plates. 

Now,  for  a  post-buckled  plate  (a>l)  we  have  restricted  ourselves  to  zero 
temperature  gradient,  i.e.,  no  thermal  loading,  so  that  the  uniform  tempera¬ 
ture  and  temperature  variation  can  be  lumped  into  the  parameter  s.  In  the 
post-buckled  state,  the  total  mean  square  displacement  is  sum  of  the  square  of 
buckled  plate  amplitude  and  mean  square  displacement  owing  to  the  acoustic 
loading  only.  As  s  becomes  large,  the  former  increases  linearly  with  s, 
whereas  the  latter  falls  off  by  1/s.  Hence,  the  total  mean  square  displace¬ 
ment  is  dominated  by  the  buckled  plate  amplitude  in  the  nigh-temperature 
limit.  It  should  be  pointed  that  the  extended  equivalent  linearization  devel¬ 
oped  in  Appendix  N  does  not  work  when  applied  to  Eq.  (7.10).  This,  together 
with  the  1/s  falloff,  is  why  the  post-buckled  plate  analysis  was  restricted  to 
zero  temperature  gradient. 
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At  this  juncture,  one  may  ask:  "How  good  is  the  random  dynamics  of  a 
single-mode  equation?"  Clearly,  we  cannot  answer  this  fully  withe Jt  the 
detailed  study  of  multi-mode  equations.  In  the  nonthermal  case,  however.  Nei 
has  shown  that  the  single-mode  dynamics  provides  an  adequate  approximation,  as 
long  as  the  external  forcing  is  weak  [11,24].  In  fact,  this  is  a  good  news  in 
that  the  thermal  effects  are  played  out  quite  poigantly  in  the  weak  external 
forcing  range.  Hence,  the  peculiarities  of  three  thermal  terms  might  have 
already  been  captured  by  the  single-mode  equation  investigated  here.  Finally, 
a  redeeming  feature  of  the  present  analysis  is  dimensionless  representation  of 
the  results,  thereby  freeing  the  validity  of  analysis  from  the  lack  of  knowl¬ 
edge  of  the  precise  temperature  dependency  of  parameters  involved  in  the 
formulation. 
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X.  Direction*  for  Further  Work 


This  report  began  with  the  derivation  of  Modal  equations  for  both  the 
simply-supported  and  clasped  plates,  under  the  immovable  edge  condition.  It 
is  this  edge  condition  that  gives  rise  to  thermal  buckling  when  the  effect  of 
large-amplitude  displacements  is  incorporated  into  the  displacement  equation. 
However,  in  an  attempt  to  exhibit  the  interplay  of  thermal  and  structural 
terms,  me  first  obtained  the  modal  equation  for  W44  ,  W  ,  W I  ,  and  K  ,  and 
then  reduced  it  further  to  the  prototype  single-mode  equation  for  quantitative 
comparisons.  Clearly,  me  notice  the  direction  of  simplification  from  multi¬ 
mode  to  the  single-mode  equation.  Hence,  the  further  mork  must  follom  the  re¬ 
versed  path  of  generalizing  the  single-mode  dynamics  to  multi-mode  equations. 

Some  specific  proposals  for  further  investigation  are  then: 

(1)  Investigate  the  dynamical  behavior  of  Eqs.  (4.4)  and  (4.11)  to  exhibit 
the  role  of  thermally  buckled  modes  on  the  modal  energy  exchange.  The  under¬ 
lying  motivation  is  that  the  post-buckled  mode  has  a  larger  amplitude  than  the 
pre-buckled.  It  is  therefore  necessary  to  investigate  the  sequence  of  modes 
undergoing  thermal  buckling  in  a  multi-mode  system. 

(2)  Validate  the  random  dynamics  of  single-mode  for  the  multi-mode  Eqs. 
(4.4)  and  (4.11).  This  should  be  carried  out  in  tmo  steps.  First,  analytical¬ 
ly  by  extending  the  equivalent  linearization  to  Eqs.  (4.4)  and  (4.11),  and 
then  numerically  by  applying  the  Monte-Carlo  computer  simulation  to  compare 
with  the  analytical  approximation.  It  is  also  suggested  that  other  random 
dynamic  techniques  such  as,  for  example,  stochastic  averaging  should  be  com¬ 
pared  with  the  equivalent  linearization  used  in  this  report. 

(3)  Investigate  t‘  &  multi-layered  composite  plates.  Although  this  is 
doable  by  the  Galerkin  type  of  analysis,  it  is  preferred  to  formulate  an 
alternative  finite  element  numerical  procedure  which  can  easily  andle  the 
complex  plate  geometry,  boundary  and  edge  conditions,  non-eymmeeric  composite 
layers,  etc. 

(4)  It  is  highly  desirous  to  provide  some  experimental  verification  for 
the  random  dynamic  analysis  presented  in  this  report.  Perhaps,  it  may  be 
necessary  to  begin  with  investigation  of  a  system  simpler  than  the  plate,  and 
a  deterministic  excitation  rather  than  the  acoustic  loading  considered  here. 
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Fig.  1  Plate  configuration. 

The  x-y  plane  Is  at  the  aid-plate  and  z-axls  Is  across  the  plate 
thickness.  The  Top  and  botton  of  plate  are  located  at  h/2  and  -h/2, 
respectively. 


Fig.  2  Coaparison  of  basis  functions  for  claaped  plate. 

(a)  Nonorthogonal  expansion  functions  of  even  aode  shape 

-  -  ♦1-2S1(x);  -  *,-2S#(x>. 

(b)  Orthonoraal  expansion  functions  of  even  aode  shape 

- f1*»/573sji(x);  — -  f#Vr2l78s#(x)  +  /37I5$t(x). 
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Configuration 


3  Typical  diitributioM  of  temperature  variation  and  gradient. 

Tq  la  dimensionless  uniform  temperature  above  the  room  temperature; 

Vv  the  *««nitude  of  temperature  variation  over  the  plate; 

V*  1#  the  *®gnitude  of  temperature  gradient  through  the  plate  thickness 


(b) 


Pig.  4  The  potential  energy  U  of  slaply-supported  plate. 

(a)  0»  (l-*)q*+|t(l~/**)(/j4+l)+2(^4+l+2p/»a)]  under  0*1  and  /**" 0.1. 

(b)  Potential  energy  surface  developing  a  symmetric  double-well. 
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■*  *  .•*=%  % 
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4 

(b) 


Pig.  8  Maximum  Mean  square  displacement  of  post-buckled  plate  (s>l). 
The  level  of  mean  square  displacement  is  raised  y  the  square  of  buckled 
plate  amplitudes,  (a)  Simply-supported  plate,  (b)  Clamped  plate. 


-50- 


3 


4 


Pig.  9  Maximum  mean  square  displacement  under  g«f-l. 

The  cross-hatched  triangle  represents  the  contribution  of  buckled  plate 
amplitude,  (a)  Simply-supported  plate,  (b)  Clamped  plate. 


Pig.  10  Total  aean  square  displaceaent. 

The  aean  square  displaceaent  Js  entirely  due  to  the  fluctuations  for  s<l, 
whereas  the  contribution  of  buckled  plate  amplitude  increases  linearly 
with  s  for  s»l.  (a)  Slaply-supoorted  plate,  (b)  Clasped  plate. 


Fig.  H  Kean  square  displacement  due  to  acoustic  loading. 

Up  to  s«l,  the  aeun  square  displacements  are  identical  to  Fig.  10. 
HoNever,  only  the  contribution  by  acoustic  loading  is  displayed  for  s>2 
(a)  Simply-supported  plate,  (b)  Clamped  plate. 
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Pig.  12  Maxi aua  aean  square  dis¬ 
placement  of  slaply- supported 
plate.  The  left  ordinate  refers  to 


the  right  ordinate  represents 
(«l£>-«^>)/<q£>  In  percent. 

(a)  TQ *0 . 2 ,  4y-«g«44  (s-0.952) 

(b)  To-0.5,  4y-4g-10.8  ( 8-0. 962) 

(a)  T  -0.8,  4-4  -2.3  (s-0.957) 


Pig.  13  Maxlaua  aean  square  dis- 
placeaent  of  claaped  plate.  The 
left  ordinate  refers  to 

-  64<q£>/9  and  -  64<q*5 

the  right  ordinate  represents 

(<qj>-<q*>)/<qj>  In  percent. 

(a)  To-0.2.  4y-4g-27  (>-0.969) 

(b)  To«0.5,  *v“V6,4  (s-0.956) 
(•)  T  “0.8,  <u»i  *1 .4  (s -0.9600 


HWM 

iVflll 


< 

a  04 


0J2  C.4  0.6  0.8 


02  0.4  0.6  3.6 


(b) 


Fig.  IS  Root-aean-aquare  of  extreae-f lber  stress  and  strain  for  claaped 
plate  under  gff-1. 

- VW0;  — 0 - V0*8*  W0; 

- T.-O.ft,  «..-6.4,  «„«0; - X -  T."0.5,  «. -«g-6.4 


Appendix  A:  Cosine  expansion  of  (|^fc)  and  expressed  by  Eg.  (3.3) 

a*«  3 

The  torn  (§]gy)  ’  Using  Bq.  (3.3)  we  have  by  the  trigonometric  identities 

(wxv)a  «  n*jT  ST  nj 
y  n*o  j-o 

tT  w|inw1jni(oos(n+i)7TX<-cos(«“i)7Tx)|(cos(n+j)ny+co8(n-j)7ry) ,  (Al) 
By  interchanging  the  order  of  summations,  the  {..}  in  Eq.  (Al)  becoaes 

{•  }  ■  £0°enj 


Here 


(A2) 


<,PnJ’5«0*,|U»W(P-“)i"<P”m)+^IpW»n,,(»-p)j"^-p'^0wBnw{B+p)j“<B+p)A<p)’ 

where  A(p)*l-4p  excludes  'the  tera  for  p*0.  (Note  that  the  factor  A(p)  can  be 
applied  to  either  the  second  or  third  tern  to  avoid  double  counting.)  Now. 
inserting  Bq.  (A2)  into  Eq.  (Al)  and  interchanging  the  order  of  suaaations, 
we  obtain  the  double  cosine  expansion 


(*_.)*=  n4  iT  if*  IL_  cospnx  cosqny, 
y  p*o  q*o  M 


(A3) 


where  6 


V*  g0°pn«I-n)n(<|-”>  *  +  £„t,pn<n.q)n<n*'»4«') 


'St  St 

The  tera  (*~)(~ ■):  In  parallel  to  Eq.  (Al),  we  have 
_ a  jr 


If  If 

xx  yy 


n4xr  r  f 

n«o  j^o 


aflf*  if*  w  w,  (-cos(»+i)irx+cos(a-i)«x)M-co8(n+j)ny+co8(n-j)Try) . 
la-o  i-o  J  J 


(A4) 


After  expressing  the  (..)  In  Eq.  (A4) 


{•  }  *  £0Snj 

where 

Cpnj“"  20*tanw(p-a)j1^  +  ^pWnnw(a-p)J^  +  ^0%nw(a+p:  J*^A(p)  * 

the  interchange  of  the  order  of  suoaations  gives 

wy  wvv  «  «4xf°  ?T  B  cospnx  cosqny, 
yy  p»o  q»o 

where 


(AS) 


(A6) 
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®pq*~  j£0Cpn(q-n)(q~n)*  +  JIqCpn<n-q)(n‘'q)*  +  £0Cpn(n^)*n+q)*A(q) 


Therefore,  by  combining  Eqs.  (A3)  and  (A6)  Me  obtain  the  final  expression 
(*xy)a_  *xx"yy  '  ^  SpqC08p"X  cosqny, 


where 


pq  pq  pq 

For  later  refi 


w31  and  ws3; 


erence ,  we  shall  enuaerate  S__  which  Involve  only  w. . ,  w.  , 

P4  X  a  X  9 


s03-2l^>-  e,o*2t<r  2»„V  sn?,- 

-  ".a"*,)- 

ea.-4f-"„»„*  O’’,,*..*  25",»»»,1-  O",,"..* 

s«-4t-9,,„’'a.1’  8..'4['9",."..'' 


Note  that  the  pattern  of  quadratic  terwa  in  the  [..]  agrees  with  Levy's 
Table  1  in  Ref  (5). 
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Appendix  B;  Evaluation  of  jy*  R4fr(x)fs(y)dxdy 

Let  us  evaluate  the  first  tens  w^F  of  K  .  In  view  of  Eq.  (3.1)  and 
(3.3),  we  find  that 


Vo  "xxFwVx>Vy>dxdy=  «*^Pxr^»P8  + 


Ehn* 


(Bl) 


where 


*r  £t  £0^£x  ^0VnPpq-,(4;-p-  6l-v+  *J+r>}<*S-q-  *q+a>- 

The  integrations  have  been  carried  out  with  the  aid  of  Eqs.  (HI)  and  (H2)  of 
Appendix  H.  By  expanding  the  factors  of  6  ,  we  see  there  are  nine  suws  in  9 


J.o^Hs-qJq  '  ^a<*atlr(q-8)q  +  ^^(q+sjq' 


where 

Zr 
p=o 

Similarly,  we  have 


q*8 


t 

(B2) 


q*o 


W  ir0*(r-P)"FM(r-p)*  -  £r«<p-r)»Wp-r>*  *  £„"<p.r),1W'm'1*  ' 


VS  wyyVr(x)lf8(y)d*dyi  n*aaPyo*wr8  + 


Ehrr4 


V 


(B3) 


where 


9  =  if*  jf°  n^/jf*  w  F  p*(4*  -  4*  +  {•  )t(jn  -  in  4,  \ 

*  n=i  q-o  '»=i  p-o  N  r"P  P"r  P+r'J'Vq  *q~®  Vs1' 


Also, 


where 


Vi  "XyFxyVx)Vy)dxdy 


Ehrr4 


(B4) 


V  £x  £<?<£>  *W  Vr>}<*S-q-  *W 

By  combining  Eqs.  (Bl),  (B3)  and  (B4),  we  obtain  the  final  expression 

Vi  *4Vx)Vy)dxdy*“  ^(b^r*  +a2py**)wrs"  *rs’ 


(B5) 


where  frs=  ft+ 29, 3  .  In  view  of  Eq.  (B2),  we  see  that  9 rs  represents  the 
nine  suns  in  Levy's  Eq.  (9)  in  Ref  [5]. 

As  an  illustration,  we  enuaerate  *„  for  w  .  w  .  w„,,  and  *„; 

*«  11  19  31  93 

fii-  Vn*  Vi.*  h,"n* 


'1.-  Vu*  VW  W  Vw 
\"ii*  k.V* 
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where 


ha-8t.r1,'u-S!F..*il'..' 

h4.-S6(Flt,F„). 
b,-100(FJ4.F„)t3«(FM*F4  J , 


h»--72(Foa+F„), 

h10"64,<,0«*r.0)-  <BT) 


Appendix  C:  Cosine  expansion  of  (|^)*and  (^!f){&~f)  expressed  by  Eq.  (3.15) 
_ y _ Ex5  dv* _ 

The  terx  (|jgfy)  •  Using  the  alternate  expression  for  Eq.  (3.15) 
"(XaV^oCo^aCn  wraaraagnSB(x)Sn(y) , 


we  have 


(wxJ*-»*2T  jf°  2Taanarfi 
*  n-o  J»o  8«n  d»j  8n 


■Co  fSt.  f!1"r*"c'i"r«aciA«<x)At(x)}An<y)*J(yl’  <cl» 

where  Aa(x)*aco3(a7rx)slnnx+sin(nnx}co3Trx.  By  trlgonoaetrlc  identities,  we 
rewrite  the  {..)  in  Eq.  (Cl) 


{**}  ”  l^isd  +  G^i#d  cos(a-i)rrx), 


(C2) 


where 


aii8d*  £ ^"cdWci2**1*1* 


r«»-«  c*j 


c-iraCdr*Cl 

By  changing  the  order  of  suaaations,  we  then  have 


to  £.H*  +  f.pa-C-P)»0  *C9*(**‘»s')4,P)5  ,C3> 

where  A(p)»(l-«£) .  Now,  substituting  Eq.  (C3)  into  Eq.  (Cl)  yields 

<■*/- Tf*  to  f-if-n  £/»^P*V»>AJ<»>}  C°*PnX'  (C4) 

where 

HP«d  ’  +  (■-*»>•<*  +  j£0<5B(a+p)8dA<p)  • 

By  repeating  the  saae  Interchange  procedures  for  the  [..]  in  Eq.  (C4),  we 
obtain 

to  £»[■  *  £,*»<»-<I>  +  ^0*pn(n»<l)4(,!)]c°*‘1"x-  (C5! 

where 

r  to  tfnKf1  poiV'FV  - 

~  ^In+#  ^j8sfn.,  )adj,,piid*n+1  H  J+1 )  ■ 

Hence,  the  final  cosine  expansion  becoaes 
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where 


(w)*-  2-  XT  jf  Z-_  cospnx  cosqtry, 
8*  p-o  q-o  M 


*  So^Mq-n)  +  £q*i>n<n-q)  +  ^^(n+q)  A(q)* 

The  tern  We  begin  with  the  representation 

_ a**. ,  .3/, 

r  ra*,, 

**  W  n-o  j-o  s-n  d-j  #n 

H.  £„{£.  Ci“i’"cd“r*°ciB,<x)Si(,‘,^s"<v’Bj<y''  ,cv 

where  D|B(x)*2aco8(*TX)cosrTX+(l+^  )Sa(x)  •  The  {..}  In  Bq.  (C7)  sunned  over  n 
and  i  gives 

£.  £„{•  }-  *£„^0a*(p-->'d  *  £p*-(»-p)*«  *  <c*> 


where 


*rs"cdarnaci2^+1)  +  (*■')* 


r-n  c-i 

+  zT  iT 

r-n+a  c« 


rs  cd  rn*ci*'"  r-n-a  c-i  r8  cd  r("‘*)  cl 

1"r."cd*r(.**)"ci(,+1>*- 


Now,  insert  Eq.  (C8)  back  into  Eq.  (C7)  to  obtain 


where 


’WW  £f-  £„  £/.n^jVdV»>BjHCM,"x’ 

*W  *  ^^(p-njsd  +  *ii(n-p)sd  +  £[0*n(n+p)sdA(p*  * 


After  interchanging  the  order  of  aunnations  for  the  [..]  in  Eq.  (C9) 
to  £,["]■  «£0t^0°P(d-J)J  *  £,°PU-q)J  *  ^PfJ-dU4"11)00'9"1'  (C10> 


where 


°pnj*  *  £„  ^J"»n"d/’p,d2(j5*1)  +  £n  ^j..,*.il!,d(j-J)np»d<J-1)* 


*  £n  £J„'*p“dU’»)T1P*d<Jfl)‘‘- 


the  final  cosine  expansion  becones 


where 


wxx"w“  “I  ^  fnoCOSrnx  cosqny, 
**  yy  8*  p-o  q-o  N 


tCll) 
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*pq  ’  +  jIq°p(J-<*)J  * $.°k j+Q)S*(q) 

Therefore,  combining  Eqs.  (C6)  and  (Cll)  gives 

^xy^8  “  "xx"yy“  n<  ^  Pna008?7”1  cosqtty, 
*  *  *  p*o  <5*0 


(C1J) 


-here  »„)/<?  ■ 

For  later  reference,  we  enumerate  P  which  involve  only  w44,  w4,,  w,4 

P4  11  X  9  *1 

and  w>#. 

P  «  3.556*?  -  3 , 975w  W  -  6.360WW  +  3.556WW  +  0.711*?, 

02  XX  XX  13  XI  9X  XX  39  19 

+3.556w4  w  -  1.272**  w  +  28.444*?-  31 .802w,4 w  +  5.689*?,, 

X J  3X  X3  33  31  91  99  93 

^20*  3  .556*^  ® '®®®**ii**i*  3,97511^^^+  3.556silWjJ+  28.444*^a 

+3.556wiaw#1-  31.802wiawJ8+  0.711**^-  1 .272wa>waa+  5.889wJ#, 

P„  — 7.111*?4  +  27.032w4w  +  27.032W4  4  w,4  ♦  14.222WW  -  18.489*?, 

22  11  11  13  11  91  IX  93  13 

-65.778W.  ,w.  +  24.487W  w  -  18.489*?+  24.487w,w-  7. 111*?., 

13  31  19  33  31  31  33  33 

P.  «- 3.558»?4  +  25.441*  w  +  6  380W44W„~  22.756w44w  +  5.689*?, 

04  11  XX  19  11  31  11  99  19 

-22.756W  w  -  10.177w4  ,W,  -  28.444*?+  203.53w,w „+  45.511*1?,, 

19  91  19  93  3X  31  99  99 

P.  *-3.556*?  +  6 . 360w  W4  +  25.441w.w-  22.756W.  4  W„-  28.444*?, 

40  11  11  is  n  si  n  ss  is 

-22.756w.,W  +  203.53W  w,  +  5.889*?. -  10. 177w,.w, ,+  45.511*?,, 

19  31  19  39  91  91  39  99 

P„ •-21 . 466w4 4  w.  +  19.20*  w„  +  19.20*?  +  19.20w.,w.  -  34.346w,,w„ 

06  XI  19  11  99  19  13  31  13  93 

-171.73*  w+  153.60W8  . 

91  99  93 

P  “-21.466WW,  4  +  19 . 20w  w„  +  19.20*  w  -  171.73WW ,.+  19.20  »?. 

•0  11  31  11  33  13  91  13  93  91 

-34 . 348«  K  +  153. 60w* 

31  33  39 

P,  *  3.556*?-  44 . 523*  w  -  23.056W4  4  *  +  29.156W..W,  -  14.2225?, 

84  XX  IX  12  11  31  IX  39  13 

+157.16wi#w81+  46.431*^*^+  17.778**^-  131 .02watw#J-  30.151w£a, 

P  «  3.556*?-  23 . Of 8w  *c  -  44.523W.  4  w,4  +  29.156ww,  ,+  17.778*?, 

42  XX  11  13  XX  91  XX  99  19 

+157.16wiSWJt-  131.02W1#WM-  14.222W^i+  46.431w#twaa~  30.151wJa, 

20.  «<*{,. 

P.  --25.60**  +  45.795w,.W  -  30.480*?,, 

90  21  91  3  99 

P_-  23.851*  w-  48.933**  w  -  21.333*?,  -  94.933**,. +  128.89W4,W,, 

26  XX  19  XX  99  19  91  13  33 

+107.81W#tWa#-  9«.427wJa, 
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P  ■  23.851*  w-  46.933i»11l»##-  94.933*f1#w#1+  107 .81»4a*a#-  21.333*8 
♦126.89*  w-  96.427*?  , 

#1  99  99 

P44»  19.081*^*^  +  19 •081*ll*#x-  34.133*^*^  +  8.533*?#-  238.93*1#*Jt 
-91.589wi#waa+  8.533*^-  91 .589*aiwaa-  27.307*^f. 

Pa»*  25  60,?a"  166  01*t,*„  +  128. 0*?,. 

P  «  25. SOW*  -  166.01*  »  +  128. 0**. , 

18  91  91  99  99 

P4f»-2.385wi;lw1#+  27.733wilwajv  2.133*J#+  104.53wi#w#J-  118.30*ia*M 
+55.335w#1w#a-  49  493*J#, 

P§4— 2,385*tl*ai  +  27.733w4J*ja+  104.53*ja*£1  +  55.335wtJW93  +  2.133*^ 
-118.30w#1n#a-  49.493*^£ , 

P_  ■- 28.80*.*..  +  25.780*'  *.„  +  25.780*.,.  *..-  23.040**  . 

M  I9g91  19  99  91  99  9S 

P4#-  137.38wi#»##+  61.440*J3,  P#4«  137.38m#1Wm+  61.440^, 


Appendix  D:  Complete  listing  of  £oa~foa 
The  enumeration  of  and  becomes  simpler  when  only  w. . ,  w...  w_4  and 

P4  P4  *1  19  #1 

***  are  retained.  As  an  illustration,  we  present  here  the  computer  listing 
of  Z _ -5„„  which  involves  a.  „  .  n_  and  only 

Oft  Oft  11  31  33 

r0 s"*oa“  32 

-48  w( 1, i)w(l,3)a(l,l)a(l,l)a(l,l)a(3,3) 

+64  w(l.l)m(1.3)a(l. 1)4(1. l)a(l.l)a(3;l) 

-64  w(l,l)w{3,l)o(l,l)a(l,l)a(3,i)a(l,l ) 

-64  w(l,l)w(3,l)a(l,l)a{l,l)a(3,3)a(l,l) 

-112  w(l ,l)w(3,3)a(l,l)a(l ,l)a(3,l)a(3,3) 

+46  w(l , l)w(3,3)a(l , l)a(l, l)a(3,3)a(3,3) 

+64  w(l,l)w(3,3)a(l,l)a(l.l)a(3,l)a(3.1) 

-48  w(l,3)w(l,3)a( l,l)a(3,3)a(l,l )a(3,l) 

+16  w(l,3)w(l,3)a(l,l)a(3.3)a(l.l)a(3.3) 

+32  w(l ,3)w(l,3)a(l,l)a(3,l)a(l,l)a(3,l) 

-112  w(l,3}w(3,l)a(i,l)a(3,3)a(3,l)a(lrl) 

+48  wil,3)w(3,l)a(l,l)a(3,3)a(3,3)a(l,l) 

+64  wCl.SjwiS.ljall.ijaiS.llaO.llaU.l! 

-180  N(l,3}N(3,3)a(l,l)a(3.3)a(3,l)a(3,l) 

+128  w(l,3)w(3,3)a(l ,l)&(3,3)a(3,3)a(3,l) 

-32  wh,3)w(3,3)a(l,l)a(3,3)af3,3)a(3,3) 

+64  w(l,3)w(3,3)a(l.l)a(3,l)e(3.1)a(3.1) 

+32  w(3, l)w(3,l)a(3,l)a(l, l)e(3, l)a(l.l) 

+160  w(3,l)w(3.1)a(3,3)a(l,l)a(3,3)a(l,l) 

-64  w(3,l)w(3,l)a(3,l)a(l.l)a(3.3)a(l,l) 

-176  w(3,l)wt3.3)a(3.i)a(l,l)a{3,l)a(3.3) 

-240  W(3,l)w(3,3)a(3,3)a(l,l)a(3.3)a(3.3) 

+416  w(3,l}w(3,3ja(3,3)a(i,lja(3,l>a{3,3j 
+64  w(3, l)w(3,3}a(3, 1 }a(l,l }a(3, l)a(3,l) 

-112  w(3,3)«(3,3)a(3,l)a(3,3)a(3,l)a(3,l) 

-272  w(3,3)w{3.3)a<3  3)a(3,3)a(3.3)a(3.1) 

+272  w(3,3jw(3,3)a(3.i)a{3,3)a(3,3)a(3,l) 

+80  w(3,3)w(3,3)a(3t3)ai3,3)a(3,3)a(3,3) 

f32  w(3,3)w(3.3)a(3,l)a{3,i)a(3,l)e(3,l).  (Dl) 
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Appendix  E:  Evaluation  of  Ksfr(x)fs(y)dxdy 

Let  us  evaluate  the  three  Integrals  In  the  order  that  they  appear  In  . 
First,  we  have 

VJ  wxxxx^r<x^s(y>dxdy* 

-n'lT  Za  Ir  f  {-4i(i*+l)C1(x)  +  ((ia+l)®+4i*  )S,  (x)  }S,  (x)dx 

■»i  l»i  j-i  “  )0  1  1  J 


The  notation  C^(x)«co8(inx)cosnx  is  used  for  the  first  equality  and  the  second 
equality  follows  froa  Eqs.  (H3)  and  (H4)  of  Appendix  B.  Using  the  notation 
£jf{I)«f(l)+f(-l)  of  Maekawa  [10],  we  have  (•  •)— S^JjIar(1+I+J)(i+I)4 . 

Hence , 

Vo  "xxxx^r  (x)f,(y)dxdy-  fir  ,B2» 

By  syaaetry,  we  can  write  at  once 

"yyyy»r<I>».<»>dId»-  ^ (E3) 

Similarly, 


Vi  *xxyyMx>Vy>dxdy 

i£4  ^1*'an{^1a»iziIjJIar(i+I+J)(1+I)*KC1anj2i2>JIa8(j+I+J)(:,+I)a}'  (E4) 


Hence,  coabinlng  Eqs.  (E2)-(E4)  yields  the  final  expresion 

/J/J  *,1’r(x)Vy)dxdy' 

Mhere  b 


(E5) 
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Appendix  F :  Evaluation  of  K4<pr(x)ya(y)dxdy 


Let  us  evaluate  the  three  teras  in  the  order  that  they  appear  in  K  . 
For  the  first  tera  we  write 

W  Nxpyy*r(x>*s<y)dxdy  -  +  V 

Here,  A,  is  evaluated  as  in  Appendix  E 

The  evaluation  of  &4  is  coapllcated  due  to  the  presence  of  Fpq 


-«r  jt  «./•  r  p  Is 

■=i  p-o  q=o  ** 

^  aacari  /(2cCc(x)-(c^+l)Sc(x))S1(x)cos(prx)dxJ- 
*^d*i  j^ndSj  {Sd(y)Sj(y)c08(<pry)dy}’ 


(PI) 


(F2) 


(F3) 


where  Cn(x)»cos(nnx)co8wx.  After  the  integration  with  the  aid  01  (H7- 

H8)  of  Appendix  H,  the  result  aay  be  put  in  the  fora 


where 


-a/8*)ir  ir  w 

n*i 


to 


^cpr"  zJSlJI(c+I)a[Ar(c+I+J-p)  +  “rtc+I+J+p)  _  ar(-c-I-J+p))  ’ 

*dqs*  -ddiZII^2~I^[as(2+I-q)  +  as(2+l+qJ  “  as(-2-I+q)] 

■<1~^d)zJZIJI(as(d+I+J-q)  +  Rs(d+I+J+q)  "  aa(-d-I-J+q)] * 


Note  that  in  the  factor  4^  picks  out  the  first  tera  d«l ,  and  only  the 
second  tern  is  retained  for  d>l.  Although  it  is  possible  to  consolidate  the 
three  terns  In  0cpr  and  Rjqg.  *»e  prefer  to  leave  then  in  their  present  fora 
for  readability. 

For  the  second  tera  we  have  by  syaaetry 

/I/;  ^pXx*r(x),»,e(y)dxdy  -  w*aaPy  +  2^  A#.  (F4) 

where 


-88- 


V  “<l/*)lT  **  if  FwP*{r"  “^cprK?  and®dqs} ’ 

*■1  n»i  p«o  q-o  '•C“i  F  ;  '•d-i  ’  J 


And,  for  the  last  tern  ws  have 

#£  ^^rtx^.tyjdxdy  -  A?. 


(F5) 


Here,  we  have 

v  <£,  £,"«&,  £"■•*«>  /*i<x)»*»(t-x)sc(,)d,} 

*{$,  J^SjNd 

using  An(x)-ncos(mrx)sinffx+sln(nnx)cosnx.  Integrating  with  the  aid  of  Eqs. 
(H9-H10)  of  Appendix  H  yields 

V  <"«*>£  £>£0  ^pQ^sAprK^.'s.Aa.}- 

where 

*ipr“sJIIJI<1+I)(ar(i+I+J-p)  '  ar(i+I+J+p)  '  ar(-i-I-J+p)]  * 

Therefore,  by  combining  Eqs,  (FI),  (F4)  and  (P8),  we  obtain  the  final 
expression 

^iVr(x)’s(y)dxdy*_  ^^x^VV'V  “  yr8-  (F®> 

where  Jr|«  A4  +  A#-  2Ap . 

Again,  as  an  illustration  t»e  enumerate  yp#  for  wlt,  wtJ|,  w#1  and  wJ#; 

v,»*  «,«„• 

J„  -  «,»„*  *,»„• 

J.i  ■  ^,1*  «.•«* 

T„  -  *,«„♦  *,«„♦  (FT) 

where 

gt -  -14 . 222(F0a+F<0  )+14 . 222(F04+F40  )+14 . 222PM-7 .  Ill  (P#4+P4a ) . 

7 . »51Poa-50 . 883P04+42 . 933Pfl#+12 . 721 ( P,0-P40 ) -27 . 032Paa+44 . B23F#4 
-23 . 851P  +23 . 056F  -19 . 081FA  >2 . 385PA. . 

at  43  44  4f 

g#-  12 . 721  (P04-P04 )  +7 . 951Fao-27 . 032FM+23 . 058Fa4-50 . 883F40+44 , 523F4t 
-19 . 081F44+42 . 933Pa(J-23 . 851Fga+2 . 385P$4 , 
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g4  -7 . 11 1  ( Fos  +Faf/ )  +45 . 51 1  ( F0 .  +F40 ) -38 . 4 ( Fog  +P#(J )  -14 . 222P,a  -29 . 156 ( P^  +Fga  ) 
+46.933(F  +F  )+34.133F  ' -27.733(F^+Pajl  ) , 

ff  -  -2 . 844F  -22 . 756FftA  76 . 8Fft. +102 . 4P„  -113 . 78( Faft -FAft ) +36 . 978Fa _ 

”8  08  04  08  08  '  80  40  '  88 

+28 . 444F. ^  +42 . 667F -51 . 2F„ .. -55 . 556F  -17 . 067F  4.267F,,.  , 

84  88  88  48  44  48 

g4»-7.111(F08+F80)+45.511(F04+P40)-38.4(Foe+Fgo)+65.778Faa-157.16(Fa4+F4a) 
+94 . 933(  Fag +F#a  ) +238 . 93F44 -104 . 53( Fgg +Fg4 ) +28 . 8Fgg , 

gL -  2 . 544F„„ +20. 353F„^ +88 . 692F„_  -91 . 589 (F„.  -F_ )+«3 . 604Fa„  -24 . 487Faa 
*7  oa  04  00  09  44  20  22 

-46 . 431F  -126  •  89FM  +166 . 01F^o  -407  .  06F^n  +131 . 02?^  +118 . 3F^m 

2*  as  28  40  42  40 

-137  * 38F_  +343 . 46F-rt  -107 . 81FM  -55 . 335F. A  -25 . 760F  +17 . 173F.0  , 

48  00  02  04  00  08 

ft  —2 . 844F  -22 . 756F„^  -76 . 8Frt  _  + 102 . 4F„„  -1 13 . 78 ( Faj>  -F^„  )  +36 . 978Faa 

8  08  04  08  08  '  80  40  88 

+28 . 444F_  +42 . 667F„ -51 . 2R,„  -35 . 556F_  -17 . 067F^  ,  4 . 267F^ , 

84  88  88  48  44  46 

g»  -  63 . 604F„„  -407 . 06F  +343 . 46F„.  +2 . 344F.„  -24 . 487F„„  +131 . 02F  -107 . 81Faa 

*9  08  04  06  80  88  84  86 

+20 . 353F_  -46 . 431F_  -55 . 335F  +68 . 692F.„  -126. 89F.a  +118 . 30F„  . 

40  48  48  60  88  84 

-25 . 760F„ -91 .589(Fa_HP  )+166.01F_  -137. 38F.a  +17 . 173F„.  . 

00  HO  44  02  04  00 

>.lo-22.756(Foa+Fao)-182.04(F04+F4O)-614.4(Foa+F#o)+818.2(Foe+Fao) 

+14 . 222Faa  +60 . 302 (F34  +F4a  )+192 . 85(Fag  +F§a  )-256(Fag  +Fga ) +54 .613F44 

+98 . 987 ( P  +F  ) -122 . 88(F,Q  +F_  )+46 . 080F„.  -30 . 720(Faa  +Foit ) .  (F8) 

44  84  t0  84  (IS  88  88 
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Appendix  G:  Evaluation  of  «(§%)*»  and  baaed  on  Eg.  (3.15) 

We  first  write  in  detail 


« 


<H>a >>*  /h|?)8dxdy-  tT  'S°  w  w  awiV|  /AjCxlA^xJdx, 

OT'  00  ”  ■*»!  n-i  r«;  *a  rnI-i  j*i  B1  rj  o  1  J 


where  A|(x)»lcos(inx)8inrrxfsln(lyrx)co8iTX.  Then,  using  Eqs.  (H3),  (E4)  and 
(H5)  of  Appendix  H,  we  obtain 

«<15>*»- Jf 

ax  8  k-i  n-i 

-  (1+l)a«r(l+2)  -  '  <01> 

Let  us  consolidate  the  {..)  in  Eq,  (Gl)  by  the  notation  Zjf fl)*f (l)+f (-1) 
of  Maekawa  [10] 


where 


«<§S>*»  - 


[G2) 


A«“"  ^,“W«>S1i»iZJIIJIar(i+W)(1+I)a' 

By  syaaetry  we  have 

<<(5y>a>>  *  "**-• 


where 


(G3) 


V" 


i £t  Ci  Ciw*nM-ianj2jFlJIa8(J+I+J5(J+I)V 
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Appendix  H:  Integrals  used  la  Qalerktn's  procedure 


/sin(ajrx)sin(bnx)cos(prrx)dx-(l/4)[dJ_p-«J_b+«®+b], 
/  cos(arrx)sin(t»rx)sin(pnx)dx«(  1/4)  [4  J_p-4  J_b-4  J+b] , 


/«in(aTrx)sin(bnx)slnanxdx*(l/8)[(2+4b)4g-4 ^_a-4^+a3 , 
/ cos  (arrx ) sin ( bnx ) cosrtx  slnnxdx- ( 1/8 )  [4^4b-4b_a+4b+a  ^ 


/cos { anx)cos ( b»rx) sinanxdx« ( 1/8)  [  (2-«b)4a- 4 a_a-<*a+a ]  , 
/sinUnxJsintbTMOcos^xdx-U/SHU^K^h,^]^] , 


/sin(arrx)8in(b7tx)cos(pnx)8ina«xoix  *(l/16)[24b_p+24a+p-24ba+p 

»b  •  b  ,b  ,b  b  ,»  b  ,»  b  i 

a->  2-p  a+2+p  a-2-p"®  a-2+p^  -a-2+p^  -a+2+p^  -a+2-p J  ’ 

1  v  w 

/C08(STTX)stn(brrx)c08(prrx)C08>TX  8iwrxdx*(  1/16) [4a+2_p+4a+2+p 

tb 


a-2-p  a-2+p  -a-2+p  -a+2+p  -a+2-p 


*-*+2-Je 


/sln(arrx)sin(brx)8in(prrx)C08sx  sinrrxdx«(l/16)[4a+2_p-4a+2+p 


-6^  «  +fib  «  -4b  «  +4b  ^  v  „  . 

a-2-p  a-2+p  -a-2+p  -a+2+p  -a+2-pJ 


/con(anx)6in(bTtx)sin(pnx)sinatrxdx«(l/16)[-24b.+24  +24b 

q  O  P  G  P  S  P 

+4b  «  -4b  o  +4b  «  -4b  o  -4b  «>  -4b  o  *4b  o  1. 

a+2-p  a+2+p  a-2-p  a-2+p  -a-2+p  -a+2+p  -a+2-pJ 


(Hi) 

(H2) 

(H3) 

(H4) 

(H5) 

(H6) 

(87) 

(H8) 

(H9) 

(H10) 


Since  the  indices  for  a  j  j  and  and  are  restricted  to  positive,  the 
teras  with  negative  indices  are  left  out  of  the  above  formulas.  For  Instance, 
a-(p+b)  *a  exc^u<*e<*  from  Eqs .  (H1-H2).  Note  also  that  Eqs.  (P  -H8)  reduce  to 
Eqs.  (H3-H4) ,  respectively,  for  p-0. 
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By  simply  ixpooirg  ajj-aj,  it  is  not  possible  to  reduce  Eq.  (PI)  of 
Appendix  F,  /J/J»»xxPyy«Pr(x)<P8(y)dx<iy.  to  Paul's  expression 

^o^o"xxFyyMxJgs*y*dxdy’  *here  frU)°cos[(r-l)jTx/a]-cos[(r+l)nx/a]  and  g8(y) 
is  defined  similarly  for  y.  The  reason  is  that  <pr  are  orthogonal  but  fr  are 
not.  To  begin  with,  let  us  reproduce  Eq.  (71)  of  Paul  [8]  in  his  notations 


^JXxFyyfr(x)*s(y)dxdy  * 

«*hbPv  r  a  a 

-  -nfir*  {-[^  l)a+(r-l)a)[(Ct+l)Wr8-  Wr(8+a)-  "r(8-a)l+ 
+(r<l)*[(Ci+l)W^r4a ja-  w(r4a){a4a)“  w(r*a)(s-a)^ 
+(r-l)a[(C1+l)W^r.t j8-  ^(r-ajjg+a)-  *(r-a )(s-a )•*} 

"  *  Wan|(*-H)a[2(n+8)aAt(n+s)-  (n+s-2)aA1  (n+s-2)~  (n+s+2)aA1  (n+s+2) 

-£(n-s)aAi  (n-s)+  (n-s-2)aAJ  (n-s-2)+  (n-s+2)aAl  (n-s+2)] 
+(m-l)a[2(n+s)aAa(n+s)-  (n+3-2}aAa(n+s-2)~  (n+3+2)aA#(n+s+2) 
-2(n-3)aAa(n-s)^  (n-s-2)aAa(n-s-2)+  (n-s+2)aAa  (n-s+2)  ]j- ,  (II) 

where 

A*  00=  C(«+r+2)F(l#+r+a)k+  C(«-r)F(1|.r)k-  C(*+r)F(l|+r)k-  C(m-r+2)F(>_r4a  )k, 
Aa(k)-  C(a+r-'2)F()i+r_ajk+  C(*-r)F(j|_r )k~  C(*+r)P(lB+r)k-  C(w-r-2)F(|J_r_a)k, 


f’2  for  k  *0 
L 1  for  k  #0’ 


/2  for  s  -1 
\1  for  s  9*1* 


We  shall  now  attempt  to  compare  Eq.  (II)  with  Eq.  (FI)  of  Appendix  F. 
First,  note  that  Eq.  (FI)  consists  of  two  terms  A#  and  A4,  given  Eqs.  (F2-F3) 
respectively.  Since  the  orthogonality  has  been  used  in  Eq.  (F2),  We  repeat 

the  derivation  of  Aa  with  the  replacement  of  /80n(y)*8(y)dy  by 

/8sin(nrry)8in(s?ry)sir?nydy  and  obtain 

n«i 

‘S.,-  ■’I-.}-  <i2 

It  is  easily  checked  that  Sw^f-.H..}  in  Eq.  (12)  is  identical  to  the  first 

(..)  term  in  Eq.  (II),  multiplied  by  it*hb?%/4ti . 

On  the  other  hand,  the  orthogonality  is  not  used  in  Eq.  (F3),  hence  A. 

■«  • 
becomes  upon  imposing  8^*4 1 
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A  *  -42f°  jT 
■*=i  n*i 

(x)-(»a+l)S||{x))S 


2T  F  qs 
P*0  t|*0 


.(x)co3(prtx)dx||/Sn(y)S8(y)co*(q7Ty)dy|, 


where  C|>(x}*co8(Rrrx)co8nx.  After  integrating  with  the  aid  of  Eqs.  (H7-M8)  of 
Appendix  H,  we  put  the  resulting  expression  in  the  fora 


A  -  -(1/8* )jf°  w__  if  F^q* 

»*i  n=i  p*o  q^o  ” 

x{‘( (a+1  )a+(a-l )a ) t< J_r+  6?(B_r)-  *S+r]+<»+l)at*S+a-r+  *Wrr 


+(a-l)  [3B_a_r+  ^-(B-sj-r)-  3a-a+r~  *-(a-2+r)ty 
x{2^n-s+  5-(n-s)_  5n+8^  “  ^n+a-s+  a-(n+a-s)“  5n+a+s^ 
~^n-a-s+  5-(n-a-s)“  3n-a+s”  tf-(n-2+s)^}' 


As  pointed  out  in  Appendix  H,  Eqs.  (H7-H8)  exclude  teras  with  negative 
indices.  Hence,  the  teras  excluded  froa  the  first  {..}  in  Eq.  (13)  are 

( (a+l)a+(a-l)a)$P(lt+r)  -  (a+l)a«?(ll+a+r),  (14) 

and 

-24Vs)  *  4-(n«..)-  I15' 

are  left  out  of  the  second  {..}  in  Bq.  (13).  Then,  inserting  Eqs.  (14-15) 
back  into  Eq.  (13)  we  have 


A  *  -(1/8*  )3f*  f  «M2r  if*  V  <f 
*  n-i  p*o  q*o  *** 

x{(a+l)  ^a+a-r+  ^-(a+a+r)+  4a-r+  a-(a-r)~  aa+r~  *-(a+r)~  *a+a-r~  a-(a+a-r)^ 
+(a-l)  (fiB_a_r+  a-(*-a+r)+  *>-r+  5E(M_r)~  tf«+r“  *-(»+r)“  4a-a-r~  *-(a+ — r)^} 
x{2tfn-s+  2*-(n-s)“  an-a+8~  tf-(n-a+s)_  ^n+a+s-  *-(n+a+s)~  23n-s"  2i-(n-s) 


+*n-a-s+  *-(n-a-8)+  an+a-s+  *-(n+2~s)}’ 


(16) 


Under  the  syaaetry  ^tp.tq^p.q  i»P*>aed  by  Paul  [6],  one  finds  that 

.){..)  in  Eq.  (16)  can  be  put  in  the  fora  of  the  second  {.,)  in  Eq. 
(11),  aultiplled  by  (T^EhVieabJDf^. 
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In  a  similar  fashion,  one  can  coapere  Bqs.  (F4-F5)  with  the  foraulas 
given  by  Bqs.  (70)  and  (72)  in  *ef.  [6]. 
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Appendix  J:  The  coefficients  at  -  aaQ 


as-(27/2)[9(l-?ia)^3|l/^+4^~*)a+l/(4^-^"1)*}+2(/»4+l+^a)]  , 

a#— 3(1  n*)fi4, 

a4-3(l-^*), 

afi  - ( 3/2 )  [{ l-/ia  )p2^9p2*4p~2+ 16/ (p+p-1  )a+l/  ( 4P+P_t  )a}+2 ( 9p4+l+10pPa )]  , 

a8-(3/2)[(l->*a)pa{s|p“3+4pa+16/(p4p~1  )z*l/(p+4p~t  )a^+2 (P4  +9+lQfjpa )j  . 

a?  *-27  ( i  - m*  )P*{p*+ 1/  ( 4P+P~X  )*} , 

a, *-27 1 1  -fi2  )p?fy “*  +1  /  (p  +4P  -1  )a} , 

a8*-3  (lV  >^{9^+64/ )a+25/  (4P+P~X  )a}, 

a10— 8(l-^)^*{v‘*+64/^’1  )a+25/(P+4P~1  )*}, 

a»i*48tJ!“it*)^a/^+^~1  >*• 

a.R-27(l-Ma)Pa{4(^a+^_a)+25/(^+Vi  )a+25/(4P+p-1  )a}, 
ai#*(3/4 ) [( 1-Ma  )P*{p*  *8  t^'a}+2(i»4+8:  *16^ )]  , 

ai4*(27/2)[(l-#ia)^a^  *36p~*+144/(ft+9p~tf+9/(4p+9p~x)*}+2(p*  +9+10(1^*)]  , 
a48— 243(l-^a)t 

ai#  •  (3/2 )  J  ( l-/ia  272/  (p*+p’a  )a+625/  (P+4P*1  )a+625/  (4 P+P"1  )a} 

+2(9P4-e^2*iPa)]  . 

ai7*(3/4;[(l-Ma)^a{8lpa+P'a}+2(8V4+1+18#^a)]  , 
a1(#— 243(l-/ia)p4. 

a18-(27/2)[(1-Ma)^a«[36/»a^"a+144/(9^*1)a+9/(V+V1)a}+2(9^4+l+10^a)]  , 
a#0 -  ( 243/4 )  [ ( 1- M* )  0**1 )  +2 (pSl+JJpp* )]  . 


For  P-1  we  find  that  l  -  -i  reduce  to  the  coefficients  C  -  C  previously 

1  SO  1  S3 


defined  by  Lee  [15] 


VV«- 


C  *a  /4.  C  *a  /4. 

9  9  3  ao 


-a#/4*-a4/4f  C8*ea/4*a#/4, 


C8-ata/4-u17/4,  c7‘*att/4,  C8-eT/4-a8/4,  C8*-a#/4«-a1S(/4, 
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Note  that  the  factor  (1/4)  Is  due  to  *  *n  contrast  to  $a«sin(airx) 

used  in  Ref.  [13]. 


Appendix  K:  Listing  of  SSP-5  under  g«1.2  end 


(SSP-5)  -  0 . 32205e+02xW( 1 , 1 )W( 1 , 1 ) 1 , 1 ) +0 . 22165e+03xW( 1 , 1 )W{3 , i )W( 3 , 1 ) 

‘ 1  -0 . 16796e+02xW ( 1 , 1 )W ( 1 . 1 ) W( 1 , 3 ) -C . 10298e+03x*f  ( 1 , 1 )*f  ( 3 . 1 )M ( 3 , 3 ) 
+0 . 15569e+03xW( 1 , 1  )W( 1 . 3)W( 1 , 3) -0 . 81000e+01x*f( 1 , 1 )W( 1 , 1  )W(3 . 1 ) 

-0.52C05e+02xW(l,l)W(1.3)W(3,3)-0.11364e-*03xK(l,3)W(3,l)l»(3,l) 

+0 . 36889e+C3xV(  1 , 3 )»( 3 , 1 )  K(  3 . 3 )  +0 . 30093e+02x*f  ( 1 , 1  )W(  1 , 3  )M  ( 3 . 1 ) 
-0 . 89215e+02xW (l,3)W(l,3)W(3,l)+0. 24040e+03xW ( 1 , 1 )W( 3  *  3)W (3 , 3) , 

{ SSP-5 )  —0 . 55987e+01xW( 1 , 1 )  W(  1 , 1 )  W(  1 . 1 )  -0 . 1 1364e+03xH( 1 , 1 ) W( 3 . 1 )W( 3 . 1 ) 

11  +o .  15569e+03xW(  1 , 2  )W(  1 , 1  )W  ( 1 , 3)  +0 . 36889e+03xtf  ( 1 , 1 )W( 3 , 1 )W (3 , 3 ) 
-0 . 26003e+02xW( 1 , 1 )W( 1 , 1 ) W( 3 , 3 ) -0 . 17843e+03x*f( 1 , 1  )W( 1 , 3)W(3 , 1 ) 

+0 . 15046e+02xW(  1 , 1 )  W(  1 , 1  )W(  3 , 1 ) +0 . 77192e+03x*f(  1 1 3 )  W(  1 , 3  )W  ( 1 , 3 ) 
+0.18377e+04xW(i ,3)W(3f 3)W(3,3)-0.65610e+03xW(l ,3)W(1 ,3)W(3,3) 
+0 . 84062e+03xK( 1 , 3 > W( 3 . 1 ) W( 3 , 1 ) , 

(SSP-5)  «-0 . 27000e+01xW( 1 , 1 )W( 1 , 1 )W( 1 , 1 )-0 . 89215e+02xt»( 1 , 1 )W( 1 , 3)W(1 , 3) 

81  +0 . 22165e+03xW( 1 , 1 )W( 1 , 1 )W( 3 . 1 ) +0 . 36889e+03xW( 1 , I )W< 1 , 3)11(3 , 3) 

-0.51491e+02xW(l,l)W(l , l)W(3,3)-0.22727e+03xW(l,l )W(1 .3)tf(3, 1) 

+0 . 15046e+02xW(  1 , 1  )W  ( 1 , 1 )  W(  1 , 3 )  +0 . 84062e+03xW ( 1 . 3  )W  ( 1 . 3  )W(3 , 1 ) 

+0.15191e+04xW(3,l)lf(3,l)W(3,l)-0.13805e+04xtl(3,l)*r(3,l)K(3.3) 

+0 . 31848e+04xW(  3 , 1 ) W(  3 , 3  )*f ( 3 , 3 ) , 

(SSP-5 )  —0 . 26003e+02xW(  1 , 1 )  W(  1 , 1  )W(  1  f  3)  +0 . 24040e+03x#(  1 . 1  )W(  1 , 1  )W(  3 . 3 ) 

**  +0 . 36889e+03xK(  1 , 1  )*f(  1 , 3  )W(3 . 1 )  -0 . 51491e+02xW(  1 , 1  )W(  1 , 1  )W(  3 , 1 ) 

-0.21870e*03xW(l,S)W(l,3)W(l,3)+O.J8377e+04xW(1.3)W(l,3)M(3,3) 

-0 . 45350e+03xW(  3 , 1  )W  ( 3 , 1  )*f(  3 . 1 ) +0 . 31848e+04xll  ( 3 . 1 )W( 3 , 1 ) W( 3 , 3 ) 
+0 . 26086e+04xW( 3 , 3 ) W(  3 , 3 )W{ 3 , 3 ) . 
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Appendix  L:  Sunnary  of  h*s  and  g's  t-nder  Ea.  (6.8) 

Under  Eq.  (6.S),  Eq.  (4.8)  reduces  to 

ht«A(l+0"*),  ha— A[i  +2/ (0* +1)3 ,  hj-AW^i/l^+DJ, 

VhiO*0’  V-h7-9A/0*.  hf-8A/(0*+l),  h#— h8-9A.  (LI) 

where  A-4vTc/2. 

On  the  other  hand,  Eq.  (F8)  of  Appendix  F  reduces  to 

gt«  14  222B(l+0~*+l/(0*+l)) , 

^=-7.951B  -12.7218/0*-  27.O32B/(0*+l) , 
g#— 12.721B  -7.951B/0*-  27 .O32B/(0*+l ) , 
g4-  7.111B(l+0~*)  -14.222B/(0*+l), 
g#-  2.844B  +113.76B/0*+  36.978B/(0*+l) . 
g#-  7.111B(l+0~*)  +65.778B/(0*+l), 
g7— 2.544B  -63.6048/0*-  24.487B/(0*+l) , 
g#-  2.844B  +113.78B/0*+  36.978B/(0*+l) , 
g8— 63.604B  -2.544B/0*-  24.487B/(0*+l) . 

g10«  22.7568(1+0“*)  +14.222B/(0*+l),  (L2) 

where  B«5yT0/16. 


i 
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Appendix  M:  Random  response  under  stationary  zero-mean  Gaussian  excitations 

For  completeness,  we  shall  present  here  the  Wiener-Khinchin  theorem,  the 
input-output  relation  for  a  damped  linear  oscillator  under  stationary  zero- 
mean  Gaussian  excitations,  and  the  equivalent  linearization  technique  applied 
to  a  Duffing  nonlinear  oscillator  under  the  same  excitations. 

Wieaer-Khinchln  theorem:  For  a  stochastic  process  x(t)  we  define  :be 
correlation  R_„(t. ,t  )»<x(t. )x(t. )>,  where  <  >  is  the  ensemble  average.  When 

••A  X  X  X  8 

the  process  is  stationary,  the  correlation  function  depends  only  on  the  time 
difference  r»t  -t  .  hence  R__(t4 ,t  )-R__(r) .  (Note  that  r  was  the  dimension- 
less  time  in  Sec.  VI.)  Moreover,  when  the  process  is  ergcdic  the  ensemble 
average  may  be  replaced  by  time  average;  however,  we  shall  use  <  >  to  denote 
both  the  ensemble  and  tine  averages.  The  correlation  Rxx(r)  is  expressed  by 
the  power  spectral  density  function  g^tf)  and  the  inverse  relation  also 
exists  through  the  Fourier  transform 

Rxx<t)-_/  gxx(f)ei2nfTdf, 

*xx(f)*  7  Rxx(T)e  l2nfTdT*  (M1) 

-ao 

Using  the  angular  frequency  «*awf  and  d^w)1  gKX(f  )/2?r,  Eq.  (Ml)  becomes 

Rxx(r)*  7 
— 

»«<"»’  5 si 

known  as  the  Wiener-Khlnchin  theorem  [8].  One  often  finds  Eq.  (M2)  defined  in 
terns  of  Gxx(<J)m^xx(c>)  •  with  the  constant  pair  (l,l/2r)  of  Eq.  (H2)  being 
replaced  by  (1/2 ,?./«)  [25].  Since  H^fr)  and  d  («)  are  even  functions,  Eq. 
(M2j  may  be  reduced  to  the  cosine  Fourier  transform  relations 

oo 

Rxx(r)«  2/  )C0SMT«Jy , 

0 

oo 

~f  R^frjcoscrdr.  (M3) 

0 

It  is  important  to  note  that  Rxx(0)*<xa>  is  the  mean  square  amplitude  (i.e., 
the  energy),  hence 

<«*>•  /  d^^)^,  (M4) 

•^0 

which  states  that  the  power  spectral  density  represents  energy 

contained  in  a  small  band  of  Au  ebout  «=±«*. 

Linear  input -output  relation:  Consider  a  damped  linear  oscillator 
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*x  +  fix  +  kx  *  f(t), 


(M5) 


where  p  is  the  damping  coefficient  and  k  the  stiffness  (Note  that  P  was  the 
aspect  ratio  in  Sec.  II).  Assume  that  f(t)  is  a  stationary  Gaussian  process 
with  the  power  spectral  density  0ff(‘>).  In  terse  of  the  frequency  response 
function  of  Eq.  (MS) 

H(«)  *  - -  ,  (H6) 

(Jc-*T)  +  ipo 

the  input-output  relation  is  [8] 

#„<«)  -  0ff(«)|H(«)|*.  (M7) 

Now,  it  is  siaplest  to  let  0ff(t>)  be  a  constant  K  over  the  range  of  a  in  which 
|H(«)j*  is  significantly  different  froa  zero.  We  then  have 


<x*>-  K  /  |H(w) |8d«. 


(M8) 


To  evaulate  the  Integral  I-J^i H(«)j*du  by  the  method  of  residues,  we  write  it 
as  I^/^f (z)dz,  where  f(z)*l/[(k-^)*+(/Jz)a).  Since  the  poles  of  f(z)  are 

z  =[N4k-0*+ip]/2,  z_«[-^4k-|9®+i^]/2,  z  »[^4k-^*-i0]/2,  z  =[-'4k-/»a-ip]/2,  (M9) 
18  #  • 

the  residues  of  the  siaple  poles  zt  and  z2  give 

I-  2ui[(z-ztf(z)iz.z  +(z-za)f(z)!z.z  ]-  (M10) 

Hence , 

<x*>*  jj|.  (Mil) 


In  standard  linear  oscillator  notations  p*- 2f«0  and  k«ua,  we  therefore  recover 

<xa>«TrS/2<’«* ,  given  by  Eq.  (5-42)  of  Lin  [8].  Since  for  saall  p,  the  I 

represents  contributions  froa  the  sharp  resonance  peaks  at  cm  t/k,  one  aay 
approximate  <x*>  by 


<x*>  * 


w0jff  (v^t) 

— W  1 


(W12) 


which  is  asyatotically  correct  as  p  *  0  (see.  Fig  5.3  in  Ref  [8]). 

Equivalent  linearization:  We  now  consider  a  damped  Duffing  oscillator  with 
the  cubic  stiffness  tern 

‘x  +  pk  *  kx  +  rx*«  f(t),  (M13) 

where  y  denoting  the  strength  of  hard  spring.  Rather  than  solving  Eq.  (M13) 
by  perturbation,  the  ala  is  to  replace  it  by  a  linear  systea  of  the  fora 
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'x  +  fix  +  kex  +  e  =  f(t). 


(M14) 


By  a  judicious  choice  of  the  oquivalont  stiffness  ke,  one  attenpts  to  capture 
the  effect  of  nonlinearity  of  Eq.  (K13)  in  a  statistical  sense,  and  the  degree 
of  failure  is  quantified  by  the  error  tera  e«(-ke+k)x+yx* ,  which  is  nothing 
but  the  difference  of  Eqs.  (N13)  and  (K14).  In  the  equivalent  linearization 
technique,  ke  is  found  by  Minimizing  the  mean  square  error,  i.e.,  d<e^>/dke»0. 
Under  the  assumption  that  x  is  Gaussian  with  zero  aean,  a  simple  expression  is 
obtained 

ke  -  k  +  37<xa>.  (M15) 

When  the  error  tera  is  suppressed  we  see  that  Eq.  (M14)  has  the  same  fora  as 
Eq.  (M5).  In  view  of  Eq.  (N12),  we  therefore  have 


<x*>  a> 


n*f  f  <v*e  * 


(K16) 


Now,  suppose  that  is  wore  or  les;.»  flat,  i.e.,  *ff (v'E) •  Then 

inserting  Eq.  (M15)  into  Eq.  (N16)  and  identifying  (JE)/flk  >8  the 
aean  square  amplitude  of  the  linear  systea,  the  aean  square  amplitude  of  Eq. 
(M13)  is  given  by  the  quadratic  equation 


3|(<j^>)a  +  <x?>  -  <>^ln>  "°- 

The  positive  root  of  Eq.  (N17) 

is  physically  lulevant  [9].  (It  is  not  necessary  to  solve  Eqs  (M15-M16) 
iteratively,  as  suggested  m  Ref  [11].) 

It  is  instructive  to  examine  the  two  limits  of  Eq.  (M18).  First,  when 
the  cubic  nonlinearity  is  weak  we  obtain  by  Taylor  expansion 


12*<xlin>  ni  .  ,  6r<xlin> 

- R -  *  1  +  R -  • 

Then.  Sq.  (M18)  becomes  <j^*e<x5jn>,  as  expected.  Second,  for  a  strong  non¬ 
linearity,  by  ignoring  +1  under  the  radical  and  -1  in  the  square  brackets,  we 
obtain 

honce  <x*>  ia  proportional  to  <xUn>lA  or  ^ff>4/*- 


(M17) 


(M18) 
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Appendix  N:  Extenalon  of  the  equivalent  linearization  technique 

Suppose  a  constant  forcing  fQ  is  added  to  the  otherwise  zero-sean  Gaussian 
excitations  for  a  damped  haraonlc  oscillator 


*x  +  fix  +  kx  *  fQ+  f(t). 


(Nl) 


where  fi  is  the  damping  coefficient  and  k  the  stiffness.  The  general  solution 
of  Eq.  (Nl)  has  the  fora 


where 


x(t)=  e"0t/2(c  co»t+C I  slriot]  +  /h(t-r ) (f  +f (r ) )dr , 
1  *  o  0 

h(t)»  Ie “^/^sinSt/S  for  t>0, 
to  for  t<0. 


(N2) 

(N3) 


and  «»|(4k?  -fP  •  (In  standard  oscillator  notations  and  k»w* ,  Bq. 

(N3)  reduces  to  h(t)«e"^wotsii£>t/C,  where  ««w()(l-<a [26J.)  The  first  term 
of  Eq.  (N2)  is  a  homogeneous  solution  reflecting  the  initial  condition  speci¬ 
fied  by  constants  and  C#,  which  eventually  dies  out  as  t-n*.  The  second 
Integral  tera  is  the  particular  solution.  Consider  the  contribution  of  fQ 


/h(t-r)f  dr  ■»  -^[-e'^t/2(fsin«t+<jcos«t)+«] .  (N4) 

o  °  ku  c 


After  a  long  tiae,  Eq.  (N4)  settles  down  to  steady  state 

«•  f. 

/h(t-r)f  dr  •*  3  , 

0  OK 

which  we  shall  denote  by  x;  i.e., 


(N5) 


(N8) 


With  this  prellainary,  we  are  now  in  a  position  to  formulate  the  input- 
output  relation.  By  splitting  x  into  the  mean  x  and  fluctuation  y 


we  rewrite  Eq.  (Nl)  as 


x  **  x  +  y. 


(H7) 


*y  +  fiV  ♦  ky  +  (kx  -  f0)  »  f(t). 


(N8) 


Because  of  Eq.  (N8),  we  find  that  Eq.  (N8)  reduces  to  Eq.  (M5)  of  Appendix  X 
after  a  long  time,  hence 


<y*>  « 


**ffU/R) 

~~n — 


(N9) 


is  the  stationary  response,  as  already  given  by  Eq.  (N12)  of  Appendix  X.  Since 
x  is  constant  and  <y  •>«;>,  we  have 
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<xa>  =  x*  +  <y*>. 


(N10) 


In  other  words,  the  total  Mean  square  aaplitude  is  sua  of  the  squared  steady 
level  and  aean  square  aaplitude  due  to  the  zero-aean  Gaussian  excitations. 

Our  goal  here  is  to  estiaate  the  response  of  a  daaped  Duffing  oscillator 
under  nonzero-aean  Gauerian  excitations,  that  is 


*x  +  fix.  +  kx  +  rx*  a  fQ+  f(t).  (Nil) 

We  again  split  x  into  the  aean  x  and  fluctuation  y,  and  recast  Bq.  (Nil)  into 
the  following  fora 


'y  +  fiy  +(k+3yl?)y+  3yx ji^+  ry*  +  kx+  rX*~  fQ"  f(t).  (N12) 

We  shall  now  replace  Eq.  (N12)  by  an  equivalent  linear  systea 

*y  +  Py  +  key  +  e  *  f(t),  (N13) 

where  the  error  tera  is 

e*(-ke+k+3yj?  )y+  3yxy*  +  7^  +  fcx+  yl?-  f Q  .  (N14) 

Following  the  equivalent  linearization  procedure,  the  ke  is  chosen  to  ainiaize 
the  aean  square  error;  i.e.,  d<ea>/dke*0.  Under  the  usual  zero-aean  Gaussian 
£.ssuaption  for  y.  we  obtain 


ke“  k  +  ayl? ♦  3y<j^>.  (N15) 

Again,  dropping  the  error  tera  in  Eq.  (N13),  we  can  write  down  the  aean  square 
of  aaplitude  y 


<y*>  as 


7T#ff  (yite) 
~&re 


(N16) 


in  analogy  to  Eq.  (N9).  Inserting  Eq.  (N15)  into  Eq.  (N16)  and  denoting  by 
<y^n>ajjt0££ (v^gJ/^k  the  aean  square  aaplitude  of  linear  equation,  the  positive 
root  of  the  quadratic  equation  is  given  by 


<y*> 


— E - 


(N17) 


in  analogy  to  Eq.  (N18)  of  Appendix  M.  Note  that  the  aean  x  in  Eq. 
be  related  to  fQ  via  Eq.  (N6)  which  now  has  the  fora 


x 


(■«7)  can 


(N18) 


8y  inserting  Eq.  (N15)  into  Eq.  (N18),  we  obtain 


* 


« 


# 


-84- 


3yx*  +  (k+3r<y*>)x  -  f#  «  0,  (N19) 

the  real  root  of  which  is  given  by 

x  =>  A  -  B,  (N20) 

where 


The  pair  of  Eqs.  (N17)  and  (N20)  can  be  solved  for  <y*>  and  x,  afterwards  the 
aean  square  of  aaplitude  x  is  obtained  by  Eq.  (N10). 

In  the  Uniting  case  of  fg*0,  Eq.  (N20)  yields  x  >"0  because  A«B.  Then, 
Eq.  (N17)  degenerates  to  Eq.  (M18)  of  Appendix  N,  which  was  derived  under  the 
assuaptlon  of  zero-aear.  Gaussian  excitations. 


» 


% 
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Appendix  0:  Noraal  stress  and  strain  components 

We  consider  here  only  the  noraal  coaponents  of  stress  and  strain  tensors. 
Let  us  begin  by  writing  the  three  teras  separately  in  the  order  that  they 
appear  in  Eq.  (2.24) 

°x  “  +  *x  +  4’ 

ay  "  a?  +  ay  +  ay’  (01) 

where  the  superscript  'a'  denotes  the  aeabrane  stress,  'b'  the  bending  stress, 
and  't'  the  theraal  stress.  By  using  the  nondlaensional  variables  of  Eqs. 

(4. 1-4. 2),  the  three  teras  in  Eq.  (01)  can  be  put  in  a  diaensionless  fora 
which  we  shall  denote  by  an  overhead  karat 


,  -0  ,  _  t 

ax  +  ax  +  V 

■A  _  -AL  A  ft 

<fy  *  <*y  +  Oy, 


(02) 


so  that  the  stress  is  now  aeasured  in  units  of  TraEha/b8;  i.e., 


a 


x 


TtW 

ba 


(03) 


Also,  by  inserting  Eq.  (03)  into  Eq.  (2.1),  we  obtain  the  expression  for 
the  strain  tensor 


.  _  n8^  7  .  „  **h*  ^ 

£  ss  •  £  £  a  —  — —  £ 

sr  .  a  » *  st  n  tt  • 


kS  X 


y  ba  r 


(04) 


Here,  the  nondlaensional  *  and  c  are  given  by 

*  y 

^  ^  ^ 
t  *  0  -  +  ■ T*T , 

x  x  ^h8 

<A  *A  A  dl  A 

«y  -  Oy  -  pox  -  ^T*T, 


(05) 


where  T*«T*  for  a  slaply-supported  plate  and  T*»T*  for  a  clasped  plate. 


A.  Expansion  for  the  stress  tensor 

Let  us  now  express  the  stress  tensor  in  series  expansions. 
Slaply-supported  plate:  Introduce  Eqs.  (3.1),  (3.3),  (3.5)  and  (3.12)  into 
Eq.  (01).  With  the  use  of  Eqs.  (4. 1-4.2),  we  obtain  the  diaensionless  stress 
coaponents  when  T’-r^h^/^+ll/lBab^l+M) 


(#*+i)r 


2  2  ^(aV+nV) 


12(l-#ia)  2(1-*")  a-i  n«i 

-  u  io(>7w),‘co*lmI  c"m 


0 
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09  09  *  g  . 

-  Z  2  <  — — “ - -^cospTTX  cosqny, 

P-o  q-o 


°x  *  7 51  z  *  (■*^*+nV)¥’B(x)V'n(y) . 

x  (1-Ma)  »-i  n-i  *”  *  n 

a1  -a*  ,  _  .(g?_tl)2»  , 

*  »  12(1-^*) 

-  (^^DT  -  «•  « 

O?  - r2  ♦ - -  Z  Z  WV  (**0*m+^ ) 

^  12(l-**a)  2(l-#*a)  »*t  n-i  *n 

-  L  co,,w 


I2(l+M)  p*0  q«0  p*0*+q* 

-  **  z  Z  / - ^ - }p* 

p-o  q-o '•(^+qa //»)*•' 


cospnx  cosqny, 


oH - *-z-  Z  Z  t»JBn(«a^*M+na)f_(x)fn(y).  (06) 

y  (l-/ia)  n-i  n-i  *“  ■  n 

where  Z-z/h  ranges  over  (1/2. -1/2). 

Clanged  plate:  Instead,  Introduce  Bqs.  (3.1),  (3.12),  (3.15)  and  (3.16)  Into 

Eq.  (01).  Then,  with  T*-nak?(/94+2Pa/3+l)/3oet?  (l^)(/»*+l)  we  obtain  the 
diaensionless  stress  components  in  a  siailar  fashion 

(/?*+20*  /3+l)T  , 

o"  - - — — — -  + - -  [/^A +M01 

x  $(l-^)tf**l)  2{l-^a)  *  8 

-  E  z  (~a7aaia)qScospnx  cos^V 
3(H#O(0*+1)  p-o  q-o 

90  OO  P  ^ 

-  Z  Z  4 — —- 22 - “>q*coapnx  cosqny, 

p-o  q-o«-(^+qa//n*J 


UVM* 


2  £  wam<Pii(*)<Pn(y)  s  w,»n9-(x)?"(y)}, 

*-»  n-i  mi  n  *-i  n-i  n  J 


't  ^t 

oj  -  oj  -  - 


t  .  _  (g*+2/r/3+nzft 
y  3(1-m*  )(/**<!> 


.  {p*+2p* /3+l)T  -  _ 

0y  ’  3(l-*i® )  *  2(1-*** )  {  PA*+A«1 


-UL+2/T/3+*)  j  j  (— )p*cospnx  cosqny 
3(l+**)(0*+l)  p-o  q-o  p 


OO  OO  r 

-  P*  E  E  /■ 


J2SL 


p=o  q=oMF^+q*/p); 


cospnx  cosqrry, 


O 


UV5) 


j^E  E  W-riP:(x)9n(y)  +  E  E  IVMxfe-fy)}.  (07) 
I**  n=i  n=i  n=i  n=i  u  J 


Under  the  temperature  variation  discussed  in  Sec.  VI  for  which  T  obey 

P«l 

Eq.  (6.8),  '*ie  sun  involving  in  o"  has  only  two  terns 
oo  oo  T  6  T  6  T 

E  E  (~a~"a^""g )qa cospnx  C08qrcy= - x-£C082»Ty+  — ~ — cos2nxcos2rry ,  (08) 

p*o  q*=o  pT+r  4(0*+l) 

and  the  corresponding  sun  in  o”  reduces  to 

oo  oo  T  $  x  i  t 

E  E  (— r-~^~r ) p? cosprrx  cosqrry* - ^-®,cos2nx+  — — cos2axcos2iry .  (09) 

p*o  q=o  P  /^+cT  4£*  4(0*+l) 

B.  Stress  tensor  involving  Wti  only 

For  the  sinplest  case  involving  only  the  W  ,  Eqs.  (06-09)  sinplify  to 

give 

Sinpiv-supported  plate: 


x  12(l-/i*)  2(1-m*)  11 

“S  ■  ,*,<*».  <») . 

x  »  18(1-**)  ’ 

;« ..  ^V’1* ,  (^i»i)b» 

y  12(  l-/t* )  2(l-**a)  11 

4..T 


— r - C082rrxco82«y>-  £-8*  cos2Try, 

(£*+!)  J  2  ii 


-  4  °){gcos2nx  +  -j~— -yCOS2nxcos2rry^-  28^tcos27rx, 


(010) 


Clanped  plate: 


> 


1 


0 


“88- 


...  * 
x  S(lV)tf*+l)  3(l-*ia)  11 

-  i£.  )/ -eo82«y  :•  — — — co&2nxco82ny| 

sci^)(^+i)  ~  <•  (ff*+i>  > 

-  f^K*  -f£cos2ny-  ff*C084ny_  _,C08gnxc£18| *2+  C032rrxcoi4^+  C08_4rrxcog2pj 

®  n\4  16  2(p+0-1)a  (£+4p~-)a  4(4p+p_1)a  J 

°x* - %■«..»!,  [ffacos2ffX*p  (y)  +  p<r  (x)cos2ny] , 

*  (l-#ia)  11  a4  1  1 


"'t 
0„=  o„ 


x  y  3(iV)(Ai)  ’ 

5t  *  i^iy 

y  3(i-#i8)(^*+d  3(i-^»a)  11 

-  Iff .  (V. 0  J~  -C08aTX  + — fff — cos2rrxc08grry\ 

3(1^)(0*+1)  4  1  16*+1)  J 


_  32^  /lcos2Trx-  1  COS47TX-  ffaco82nxco82.iy  ,/»*Cos2Ttxco84Tty «gf coagTrxcos^ 
9  14t4  16  2(£+0-t)a  Aftf+4fl_1^a  Utf+«-M* 


4(ff+4/J  M* 


(4W  x)‘ 


y  d^) 


~Wila11  [^MCosJhrxq^  (y)+  <pt  (x)cos2*y] 


(Oil) 


Note  that  for  4V“0  the  o*  and  o*  given  by  Eq.  (010)  agree  with  Eq.  (34)  of 
Wilcox  ft  Cleaner  [16]  when  kx»ky*»  .Is  imposed  in  their  equation. 

For  computation,  however,  it  is  more  convenient  to  regroup  Eqs,  (010)  and 
(Oil)  in  powers  of  WljL 
Siaply-supported  plate: 


0^+1  )Tr 


■ - ~fl  +(  1-*m)  (— tt)  f~COS2rry  +  — | - COS2nxcos2ny]  +Z4_siimxsimry| 

12(l-/ia)L  1  (0*+l)  J  8  J 

♦  a i nit xs i tm yW  +  If  (ff2^)  _  0*cos2rry\w“  . 

(.-/»*)  11  ^-(1-M8)  i  11 


“ — - — - — r^fl  *(1-m) (~ri)f_C082ftx  +  — cos2*xco82jry]  +Z4 -SinRXsinryl 

y  12(  l-**a  )J-  4  *•  (/S*+l)  J  *  J 


(1-M*) 

Clamped  plate: 


^p^i-sinrrxsinnyW  +  -  cosZnxlw*  . 

-u*l  11  (1-M  )  ■>  11 


(012) 


(?*+2ff*/3+l)T0 

3(1-Pa)(^*+1) 
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x|l  +  (l-/i)(-j)^C0s2rry  +  — -*c0827Txco82rryj  +2^^3ina7TX8insrry^ 

- -[^cosZnxsit^Tiy  +  p8il^7TXCOs2»ty3W  + 

3(1-P8)  11  3(1-m*)  11 

P#£cos2ny-  ftcos4ny-  -£222S2£2S2El  +  cos2ffxco.34ny  f  cos4Ttxco32rry>  ^ 

2 (fi+P~1)9  4(40+/T1)8  J  11 ' 


(£4+2P*/3+l)T0 

y  8(lV)(^*l) 

»[l  *  (l-#i 

- ^r-t^MCOSZrrxsil^ny  +  sir^nxcos2iry3W  + 

3 ( 1 “M8 )  11  3(1-m8)  11 


)(— J)^-C032nx  +  -Jl* — COS&txCOs2»tyJ  +Z6K8in8rrX8lr?rry^ 


(013) 

C.  Strain  tensor  involving  Wj4  only 

In  view  of  Eqs.  (6. 4-6. 5),  the  substitution  of  Eqs.  (012-013)  into  Eq. 

(05)  yields  the  dimensionless  strain  tensor 
Siwply-supported  plate: 


*x-  iC’i^!)gvTo{fv~  ^(-(cosany-^cosarx)  +  -|^j^^cos2»rxcos2rryJ J- 
+  2Z/J8sin7rxsinTry¥it+  ||/?8-(/?8cos2ny-pcos2TTX)|w^4  , 


cy=  ~}f(I^!)gvTo{fv~  -3(-(cos27ry-^cos27rx)  +  -^^jcos2»rxcos2>TyJ J. 

+  2ZsinnxsimryWi4+  ■||l  -(cos2rry-  ji^cosETrx)^*^  , 

Clawped  plate: 


(014) 


*x" 


ItL—iiHl )&  T  / f  -  4f-(cos2ny-/icos2nx)  +  ^cos27Txeos2rcyl\ 

3(  1+M ) (^8+l )  v  0 (•  v  4<-  (^+1)  JJ 

-  ■^|^cos2«xsinsnyWii+  t 

-  ?pf4(P;,cos27Ty-/icos27TX)-  ■T7c(/»8cos4Try-pco34TTX)-  SSS^SSSSS^DL { l-jtfi2 ) 
y  1  lb  2(^-1)8 

,  cos2Trxcos47ry^1  vS2 .  ,  cos4TTXCos27Ty.  1  -a 

)8  *  (40+£“1)8  ^  J  n’ 

£v*=  LfL±M-J£±l ) g  T  / f  -  4f-(cos2ttX-MCOs2iry)  +  -^^^cosZirxcosZnyll 
y  3(1+ji)(0*+1)  v  o\  v  41  ^  (^a+1)  M/ 
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-  -^Psi^nxcosgrryW  .  4  -fwf 

w  11  Oil 

-  -|2|-|(cos2<Tx-/*02cos2rry)-  -jg(cos4rx-/i^acos4ny)- 

x  cosanxcos^V/P3  .  +  co34ffXcos2nyfga 

-i'2  4  (4/9+0  _1)2  4'->  11 


cos2nxcos2n-y 
2(0 +5  _1  )3 


(0+40  1  ) 


Here,  fv=sin*Trxsin*rry  should  be  inserted  in  Eqs.  (014-015). 


b 


* 


iP*-v) 

(C15) 


-91- 


U  S  Qov«fnm«nt  Printing  Ofl>c*:1992— *<8-127/62507 


